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Introduction 


In this unit we discuss the use of graph theory to solve certain physical 
network problems. Most of the discussion is in terms of electrical networks, 
for which the method we give can be used to find all the currents and 
voltages in the network. Electric currents and voltages are examples of 
through and across variables, and they satisfy the well-known Kirchhoff's 
current and voltage laws. The method we describe is very general, and can 
be applied to any physical network for which the through and across 
variables satisfy analogues of Kirchhoff’s laws. Examples include 
hydraulic networks, mechanical systems, pneumatic networks and 
thermal networks. 


A physical network can be represented by a graph in which the edges 
represent the components of the network and the vertices correspond to the 
terminals at which the components are connected together. Associated 
with each edge are the edge-variables — the through and across variable 
(currents and voltages in the case of an electrical network) which satisfy 
Kirchhoff-type laws. Since we need to be able to determine the directions 
of the through and across variables, we assign a reference direction to 
each edge of the graph. We call a graph with reference directions an 
oriented graph. 


If the physical behaviour of the components can be modelled by equations 
(called the component equations) then these equations, together with the 
equations resulting from applying Kirchhoff-type laws, can be solved to 
find the values of all the through and across variables for the network. 
For networks of small size, these equations can be solved by hand, but for 
larger networks it is more appropriate to use a computer. We show how 
the equations can be formulated as matrix equations; this procedure is 
appropriate for solving network problems on a computer. The required 
number of equations resulting from Kirchhoff’s laws are obtained directly 
from the incidence matrix of the oriented graph. This means that the only 
information which needs to be fed into the computer is the incidence 
matrix together with the component equations. This is particularly 
convenient, because the incidence matrix can be written down directly from 
the oriented graph. 


In Section 1, Modelling physical networks, we explain what is meant by 
through and across variables, giving examples for different physical 
systems, and show how a physical system can be represented by an 
oriented graph. 


In Section 2, Electrical networks: matrix equations, we show how to use a 
spanning tree of an oriented graph to obtain the maximum possible number 
of equations containing no redundant information by applying Kirchhoff’s 
laws to the fundamental cutsets and cycles with respect to the chosen 
spanning tree. We then show how these equations in matrix form can be 
obtained directly from the incidence matrix of the oriented graph. The 
solution of these equations for electrical networks is discussed in Section 3, 
Electrical networks: solving the network equations; in particular, we discuss 
the method of Gaussian elimination. 


1 Modelling physical networks 


A physical network can be considered to be a number of individual 
components connected together. In the case of an electrical network, these 
components might be resistors, capacitors, transistors, transformers and so 
on. We first look at the representation of individual components, and the 
mathematical modelling of their behaviour, before going on to show how 
the whole connected network of components can be represented by an 
oriented graph. 


1.1 Components and through and across variables 


Components 


We assume that the physical systems we wish to analyse can be 
considered to consist of a number of interconnected components. 


The point at which one component is connected to another is called a 
terminal. The following diagrams represent components with two, three 
and four terminals. 


A 


— 
2 < D c 

A 2-terminal component might be a spring, a resistor, a length of wire or 
Piping, or any other physical object which fits into a system via two 
connections. For example, a spring functions when both ends are 
appropriately connected so that it can be compressed or stretched. 
Similarly, a resistor functions when it is connected so that current may 
flow through it. 


We use the following diagrams to denote electrical 2-terminal 
components: 


resistor. —{———} voltage source 1O 
capacitor == + current source ‘8 


inductor, 7 


For example, diagram (a) below is the circuit diagram of an electrical 
circuit containing a voltage source connected to a resistor and a capacitor, 
and diagram (b) is the circuit diagram of a current source connected to a 
resistor and two inductors. 


(a) (b) 


A 3-terminal component might be a transistor. We discuss transistors in 
detail later in this section. 


A 4-terminal component might be a transformer. 


Ao- oC 
— >-—°D 


transformer 


The much simplified diagram above shows a transformer in which two 
coils of wire, called the primary and secondary, are wound on a magnetic 
core so that they are in close proximity to each other but are not in 
electrical contact. The letters A and B denote the terminals of the primary 
winding, and C and D denote the terminals of the secondary winding. 


We use a schematic diagram to show how the components of a particular 
system are connected together. For example, the following schematic 
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Note that for some applications one 
of the input terminals, A or B, may be 
connected to one of the output 
terminals, C or D, so that the 
transformer operates as a 3-terminal 
component, See, for example, 
Exercise 1.2(c). 


diagram represents a system comprising five 2-terminal components, two 
3-terminal components and one 4-terminal component: 


Through and across variables 


In Networks 1 we discussed flows in networks where the flow might, for 
example, be the flow of traffic in a road network. A flow in a network has 
the property that the amount of flow of a particular commodity into any 
vertex (other than a source or sink) is equal to the amount of flow out of it, 
A variable with this property is classified as a through variable because it 
represents an actual flow through part of a system. 


Many types of network also involve a different kind of variable called an 
across variable. To make clear the distinction between these two types of 
variable, let us look at a specific example, an electrical network. 


In an electrical network we are concerned with the amount of current 
which passes through a component in the network. This current is a flow, 
or a through variable, which is measured by a meter in series with the 
component as in figure (a) below. But there is also an across variable which 
is the potential difference (or voltage) between (or across) any two points of 
the network. The potential difference can be measured by connecting a 
voltmeter across the two points under consideration, in parallel with the 
component, as in figure (b) below. 
flow of current potential difference 
from A and B between A and B 
4 Q B 4 Bg 
ammeter in series 
voltmeter in parallel 
through variable across variable 


(a) (b) 


An analogous example is that of a hydraulic network, which consists of a 
network of pipes carrying water or some other fluid. Here the through 
variable is the amount of water flowing through any part of the network 
and the across variable is the pressure difference between any two points 
of the network. 


Consider now a mechanical system such as a car suspension. Here we have 
a mass (the car) supported by a spring and a shock absorber connected onto 
the wheel. As far as vertical motion is concerned, the wheel is more or less 
‘fixed’ to the road. A schematic diagram is given in the margin. 


An engineer designing the suspension system is concerned both with the 
forces in the spring and shock absorber and with the velocity of the 
vertical motion of the car. The force in the spring can be measured by 
placing a load cell (a sophisticated version of bathroom scales) between 
the car and the spring, or the spring and the road; the force is then carried 
by the load cell as well as by the spring (just as the bathroom scales bear 
your weight when you are standing on them). The force is thus measured 
by placing a meter in series with the component, and so is a through 
variable. The vertical velocity of the car can be measured, using some 
electronic timing mechanism, by measuring the change of the distance 
between a point on the car and the road in a given time: the measurement 


mass 


spring 


n 


ad 


shock 
‘absorber 


here is across the spring, from the car at one end to the ground at the other. 
Hence the velocity is an across variable. 


In each of the above examples, there are two variables: the across 
variable, representing the difference of a quantity between one terminal of 
a component and another, and the through variable, representing the 
quantity flowing through the component. In the following table we list 
some examples of through and across variables for different types of 
physical system. 


Table 1.1 Some examples of through and across variables 


type of through across 
system example variable variable 
resistor 
current, potential 
electrical va pa current difference 
‘voltage’ 
potential difference eo, 
‘voltage’ Ninos 
| 
spring 
aoe extension 
mechanical j | force =e 
x 2 ba as 
extension 
valve 
Bee ~ , fluid 
| flow rate fluid pressure 
hydraulic iH i Rk difference 
1 2 flow rate 
P,-P, 
pressure difference 
insulating 
layer 
heat energy peiperanca 
heat ene: differen: 
thermal bai BEp flowraie’ flow rate er sd 
1 
i eae 


temperature difference 


Equations relating through and across variables 


For electrical networks, the currents and voltages satisfy the following 
two laws. 


Kirchhoff’s current law 

The current flow into any vertex of an electrical network is equal to the 
current flow out of it. 

Kirchhoff’s voltage law 


The algebraic sum of the potential differences across all the components 


around any circuit (cycle) in an electrical network is zero. 


Analogues of both these laws can be given for many types of network 
containing through and across variables. For example, the analogue of 
Kirchhoff’s current law for a network carrying a flow is simply the flow 
conservation law: what flows into a vertex of a network (other than a 
source or sink) must flow out of it. The analogue of Kirchhoff’s voltage law 
is the potential difference conservation law: the sum of the potential 
differences around any circuit of a network is zero. 


G. R. Kirchhoff (1824-1887) 

Gustav Kirchhoff formulated his laws 
of electrical networks while a student 
at the University of Kénigsberg. Later 
he became professor of physics at 
Breslau, Heidelberg and Berlin. 


Example 1.1 

Consider the following diagram of a hydraulic network in which the 
number next to each vertex represents the pressure of the water or other 
fluid at that vertex. 


The pressure difference (or potential difference) across AC is equal to: 
(the pressure at C) — (the pressure at A) = Po - P, =70-50=20. 
Similarly, 
the pressure difference across CD is Pp ~ Po = 90-70 = 20; 
the pressure difference across DA is P, — Pp = 50-90 = 40. 


Adding these numbers together, we see that the sum of the pressure 
differences around the circuit ACDA is zero. In fact, the sum of the 
pressure differences around any circuit is zero — this is the analogue of 
Kirchhoff’s voltage law for hydraulic networks. a 


In the next section we show how each component of a physical system — 
and hence the whole network — can be represented by an ‘oriented’ graph. 
With each edge of the graph we can associate a through variable and an 
across variable. These variables are related by the ‘component equations’. 


The component equations relate the through and across variables for a 
single component. For example, if the component is an electrical resistor, 
the appropriate component equation is Ohm’s law: 


voltage 


it = 4 constant (the resistance). 


For a spring, the corresponding component equation is Hooke’s law: 
tension 


extension = 2 constant (the stiffness). 


In this unit we consider only systems which obey Kirchhoff’s laws or their 
analogues. We can use these laws to derive two further sets of equations, 
one set relating the through variables in the system (the vertex law 
equations) and the other relating the across variables (the cycle law 
equations). 


Given these three sets of equations — the component equations, the vertex 
law equations, and the cycle law equations — we can set about solving 
them. That is, if we know the values of some of the edge variables of our 
graph (say, the voltages at particular points in an electrical network), 
then we can work out the values of all the other edge variables — the 
voltages and currents in this case. The aim of this unit is to describe how 
we tackle this type of problem. 


Generally, we shall use electrical networks as examples to illustrate the 
solution techniques. You should, however, always bear in mind that the 
methods we describe are completely general, in that they can be used to 
work out the values of edge variables in any physical system that obeys 
Kirchhoff’s laws. When we have worked out values for the edge 
variables, we may then interpret them as currents and voltages, flows and 
pressures, forces and velocities, and so on, depending on the type of system 
in hand. 


1.2 Graphical representation of components 7, 


Graphical representation of 2-terminal components 


In this section we show how 2-terminal components can be represented 
graphically. Later in this section we show that a component with more 
than two terminals (a multi-terminal component) can be represented 
graphically by treating it as a number of 2-terminal components. 


A 2-terminal component can be represented by two vertices joined by a 
single edge. This is illustrated in the following diagram which shows a 
resistor and its representation as the edge of a graph: 


v(t), it) 
6} ¢ «13 
component edge 


Associated with the edge of the graph are two edge variables — one 
through variable and one across variable. In the case of the resistor shown 
above, the across variable is the voltage v, (t) across the resistor, and the 
through variable is the current i,,(t) flowing through the resistor. We use 
the subscript n when we wish to indicate that the resistor is the nth 
component of a particular system. To avoid cluttered diagrams, we label 
the corresponding edge simply by the letter n. 


We also need a reference direction for the edge variables. We get this by 
assigning an arbitrary direction to each edge, and we indicate it by 
drawing an arrow alongside the edge. Currents which actually flow 
through the component corresponding to the edge in the direction of the 
arrow are regarded as positive, and currents in the opposite direction are 
regarded as negative. Similarly, voltage differences between the 
terminals corresponding to two vertices are always taken as the voltage at 
the tail of the arrow minus the voltage at the head. Thus if the voltage 
at the tail is greater than the voltage at the head, the voltage difference 
is positive: if the voltage at the tail is less, then the voltage difference 
across the edge is negative. It is important to note that the reference 
direction is purely arbitrary and does not necessarily represent the actual 
direction of the current flow or voltage difference. 


Example 1.2 


Consider a system consisting of two 2-terminal components — a battery and 
a resistor. The following figure shows the circuit diagram and its 
graphical representation. 


a 


B 
circuit diagram graphical representation 


The circuit diagram shows the actual direction of flow of the current I and 
the actual direction of the voltage V across the battery. In the graph, the 
edge labelled 1 corresponds to the battery, and the edge labelled 2 
corresponds to the resistor. Since both components are connected together 
at terminals A and B in the circuit diagram, the corresponding edges are 
connected to the same vertices in the graph. 

The edge variables associated with edge 1 are the current i, flowing 
through the battery, and the voltage v, across the battery. Similarly, the 
variables associated with edge 2 are the current i, flowing through the 
resistor, and the voltage v, across the resistor. 
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We have introduced the letter ¢ to 
tefer to time, since the currents and 
voltages are often time-dependent. 
However, we usually omit t and write 
just v,, andi, 


Note that our convention of using the 
reference arrow to specify the 
direction of a voltage is the opposite 
to that normally used in electrical 
circuit theory: an arrow ona circuit 
diagram indicates a positive voltage if 
the voltage at the head is greater than 
the voltage at the tail. Here we use a 
different convention because we are 
using a single reference arrow to 
indicate the direction of both current 
and voltage. This usage is illustrated 
in Example 1,2. 


Consider edge 1. The reference direction is the same as the direction of 
flow of the current through the battery, so we have 

i, =I. 
In the circuit diagram, the arrow by the battery points towards terminal 
A, indicating that A is the positive terminal of the battery, Therefore the 
potential at terminal A is V volts higher than the potential at terminal 
B. In the graph representation, the arrow associated with edge 1 points 
from B towards A. The voltage v, is equal to the voltage at the tail of the 
arrow (B) minus the voltage at the head (A) — that is, 

vp, =-V. 


Now consider edge 2. The reference direction is again the same as the 
direction of the flow of the current through the resistor, so we have 

n=l 
Terminal A of the resistor is connected to the positive battery terminal 
and so is at a higher potential than terminal B. Thus the voltage at the 
tail of the reference arrow on the graph is greater than the voltage at the 
head, so 

v= V. 
We have expressed all four of the edge variables in terms of the two 
quantities given on the original circuit diagram. a 


We have now given examples of the graphical representation of a single 
2-terminal component, and of two such components connected together. We 
can extend this method of representation to more complicated networks. In 
fact, any system consisting entirely of 2-terminal components can be 
represented by a graph whose edges correspond to the components and whose 
vertices correspond to the terminals, 


Example 1.3 

The following electrical network contains nine components — an 
independent voltage source, two resistors, two capacitors and four 
inductors. Replacing each of these nine components by an edge, and 
assigning an arbitrary direction to each of these edges, we obtain the 
corresponding graph: 


electrical network graphical representation 


You may be wondering why we have not drawn the direction arrows on the 
edges (rather than beside them), referred to the edges as ‘arcs’, and called 
the diagram a ‘digraph’, rather than a graph. The reason is that for much 
of our discussion the directions are irrelevant: for example, we shall need 
to refer to cycles such as BCEB or BCDEB, and these can be considered as 
cycles only when we ignore the directions of the arrows. A diagram of this 
kind, which is something between a graph and a digraph, is called an 
oriented graph. 


Oriented graphs can also be used to represent systems other than electrical 
systems. For example, the above oriented graph also corresponds to the 
following hypothetical hydraulic power system, in which the terminals 
A, B, C, D, E and F are pipe junctions, and each segment of pipe contains a 
valve. The ‘voltages’ v,,..., ¥ are across variables representing the 
pressure drops across the valves, and the ‘currents’ i,, ..., iy are through 
variables representing the flows through the corresponding valves. 


Note that, although the reference 


arrow on the graph points in the same 


direction as the arrow on the circuit 
diagram, the expression for v, has a 


minus sign — this is because the sign 


conventions are different in the two 
cases. 


ll 


3 6 
1 2 5 7 
pane valve 
D 
Pe En Bmp! 
hydraulic power system oriented graph 


Problem 1.1 


Label the components and the terminals of the following electrical 
network, and draw an oriented graph which represents it. 


Component equations for 2-terminal components 


We now turn to the equations relating the two edge variables for each 
component — the component equations. Suppose that we measure the 
voltage across a 2-terminal component such as a resistor for different 
values of current flowing through it. We can then try to model the 
relationship between the voltage and current by an equation. Such an 
equation is called a component equation since it is characteristic of the 
component itself, as opposed to Kirchhoff’s laws, which are a general 
property of any electrical network. 


We start by looking at four types of electrical component — resistors, 
capacitors, inductors and independent sources. 


A resistor is modelled by a 2-terminal component whose edge variables v 
and i satisfy a relationship of the form 


v=Ri, 
where R is a constant, usually called the resistance. 


This relationship is known as Ohm’s law. It is sometimes written in the 
form i = Gu, where G (= 1/R) is the conductance of the resistor. 


A capacitor is modelled by a 2-terminal component whose edge variables 
v and i satisfy a relationship of the form 


where C is a constant, usually called the capacitance. 


An inductor is modelled by a 2-terminal component whose edge variables 
v and i satisfy a relationship of the form 


di 
v=Liy, 
where L is a constant, usually called the inductance. 


Note that the component equations for a resistor, a capacitor and an 
inductor are all linear equations. Note also that R, C and L are constants — 
we are assuming that the behaviour of the components is time-invariant. 
The concepts of linearity and time-invariance are discussed briefly at the 
end of this section. 
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To model components such as batteries and electric generators, we use the 
concept of an ideal independent source. This is an idealized 2-terminal 
component which keeps the value of one of the edge variables at a 
prescribed level independent of the value of the other edge variable. An 
ideal independent voltage source, for example, maintains a prescribed 
voltage across its terminals, irrespective of the size of the current flowing 
between them. Under certain conditions, a battery can be modelled in 
terms of such a source. Similarly, an ideal independent current source 
maintains a prescribed current irrespective of the value of the voltage 
across its terminals. 


The prescribed level of an independent source may not necessarily be 
constant; for example, for an independent voltage source, the prescribed 
voltage may be a sinusoidal voltage of the form v, sin wt, where v, and @ 
are constants. 


Some 2-terminal components of other types of system such as mechanical 
or hydraulic systems can be considered as analogues of these electrical 
components. The electrical analogues of two components of a mechanical 
system are considered in the following example and problem. 


Example 1.4 
Suppose that a force f accelerates a body of mass m. Then, according to 
Newton’s second law of motion, the acceleration a is given by the equation 
f=ma 
or 4 
uM 
famip: 


where u is the velocity of the mass. In this case, f is the through variable 
and u is the across variable. This equation is of the same form as the 
equation 

¥ dv 

i=Ci 
for a capacitor. We may therefore consider the force f to be analogous to 
the current i (also a through variable), the velocity u to be analogous to the 
voltage v (also an across variable), and the mass m to be analogous to the 
capacitance C. The following table shows the analogous quantities for the 
two components. 


mechanical electrical 
component mass capacitor 
through variable fy i 
across variable u v 
d . dl 
component equation f=m eB i=C a a 


Problem 1.2 


For a spring which obeys Hooke’s law, the extension x of the spring 
produced by a force fis given by the equation 


fakx, 
where k is a constant (the stiffness of the spring). If the rate of change of 
the extension of the spring (that is, the velocity u = 3) is considered to be 


analogous to a voltage, which electrical component is analogous to the 
spring? Construct a table of analogous quantities similar to the above 
table. 
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For the remainder of this unit, all the examples we consider are electrical. 
We can obtain corresponding results for other physical systems by 
considering the appropriate electrical analogues. 


Multi-terminal components 


A multi-terminal component is a component which has more than two 

terminals, and which cannot therefore be represented by a single edge 

joining two vertices. Before showing how such a component can be 

represented graphically, we consider two important types of electrical _ It is not necessary to be familiar with 
multi-terminal component — transformers and transistors. ) the technical details given below. 


Transformers 


A transformer is a device used extensively for changing the voltage at 
which electrical energy is supplied, A schematic diagram of a transformer 
consisting of two coils wound on an iron core is shown below. 


mturns ironcore — n turns ideal transformer 


Under certain circumstances a transformer can be modelled by an ideal 
transformer which satisfies ‘ideal’ conditions such as having zero- 
resistance windings and no magnetic field losses. 


A circuit diagram of an ideal transformer is shown above. Notice that one 
of the input terminals is connected to one of the output terminals; a 
transformer connected in this way can be regarded as a 3-terminal 
component. The equations relating the voltages and currents shown on the 
diagram, as derived from electromagnetic theory, are 
% Ae Note that v,,i,0, and i, are usually 


UV, Ny bony’ functions of time. 


where n, and n, are the numbers of turns in the two coils. 


Transistors 


A transistor has three external connections called the emitter, the base and 
the collector. Under certain circumstances, its behaviour can be modelled by 
the circuit shown in the following diagram. 


collector 
base 
emitter 
emitter 
(a) pnp transistor (b) controlled source representation. 


8mis the transfer conductance 
or ‘mutual conductance’ 
of the transistor * 


Diagram (a) shows the symbol commonly used to represent a pnp type of 


transistor. Diagram (b) shows the controlled source representation. When A controlled source is an independent 
a small voltage v,, is applied between the emitter and the base, under source whose output level is not 


certain circumstances a greatly amplified voltage will appear between Constant but is controlled by the value 
the emitter and the collector; in other words, the transistor acts as an 4 particular voltage or current. 
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amplifier. This is represented in diagram (b) by a current source which is 
controlled by the value of the input voltage v,,. The current flowing in the 
source is equal to g,, Up, where g,, is a constant for particular operating 
conditions of the transistor. Notice that resistors have been added to the 
controlled source representation to model the input resistance r, and the 
output resistance r, of the transistor. We can use this representation to 
obtain relations between the voltages and currents. Since the input circuit 
consists of just a resistance r,, the input current i, is given by 

Poe . 

" 

The output current i, is equal to the source current g,, Upe minus the current 
flowing through the resistor r,, that is, 


i,= 


v, 
ts = 
5o= 8m Ube 13! 


where v,, is the output voltage between the collector and the emitter. 


Graphical representation of multi-terminal components 


We have seen how a 2-terminal component can be represented graphically 
by a single edge. We now show that a multi-terminal component can be 
regarded as a number of 2-terminal components, and can therefore be 
represented by a graph with a number of edges. To do this, we need to 
make use of Kirchhoff’s current and voltage laws which we stated earlier, 
and repeat here for convenience. 


Kirchhoff’s current law 

The current flow into any vertex of an electrical network is equal to the 
current flow out of it. 

Kirchhoff’s voltage law 


The algebraic sum of the potential differences across all the components 
around any circuit (cycle) in an electrical network is zero. 


These two laws can be formulated in graph-theoretical terms as follows. 


l 


Vertex law 


The algebraic sum of the through variables associated with the edges 
at each vertex of an oriented graph is zero. 


Here ‘the algebraic sum’ means that the orientation of the edges must be 
taken into account, so if edges oriented into a vertex are counted as positive, 
then edges oriented out of a vertex are counted as negative. 


Cycle law 
The algebraic sum of the across variables associated with the edges in 


any cycle of an oriented graph is zero. 


Here ‘the algebraic sum’ means that we choose a particular direction 
around a cycle (say clockwise) and count edges oriented in the same 
direction as positive and edges oriented in the opposite direction as negative. 


We know that a 2-terminal component can be represented by a single edge. 
Associated with the edge are two variables: the voltage across the 
terminals of the component (across variable) and the current flowing 


i) tig + y-ig=0 


0s, % 


D3 
01-0) + 05+ 4-05 =0 


15 


through the component (through variable). We now consider the 
variables associated with a 3-terminal component. 


We can measure the current flowing out of each of the three terminals A, B 
and C; the currents are shown on the diagram as i,, ig and ic. Similarly, 
we can measure the voltage across each pair of terminals; the voltages are 
shown as vg for terminals A and B, vg¢ for terminals B and C, and vc, for 
terminals C and A. 


At first sight it would seem that we need six variables to characterize 
a 3-terminal component — three voltages and three currents. However, 
these six variables are not all independent since we can obtain 
relationships between them by applying Kirchhoff’s two laws, as 
follows. 


Kirchhoff's current law states that the total current flowing into a vertex 
is equal to the total current flowing out of the vertex. It can be deduced 
that the total current flowing into the region bounded by broken lines in 
the figure below is equal to the total current flowing out of that region. 


Applying Kirchhoff’s current law in this way to the currents flowing 
across the boundary in the figure, we obtain the equation 


in tig tic=0 


ip =-ip-ic. 
A consequence of this last equation is that if we measure ig and ic then we 
do not need to measure i4; we can calculate it from the equation. So only 
two of the three current variables are independent. 


We can apply Kirchhoff’s voltage law to the cycle ABCA. We thus obtain 
the equation 


Vap + Upc + Vc, =0 


Yc ™=—VaB- ca: 
Thus if we measure v4, and vc,, then we do not need to measure vgc; we 


can calculate it from the equation. So only two of the three voltage 
variables are independent, 

Thus, of the six current and voltage variables, only four are independent — 
two current variables and two voltage variables. It follows that, to 
describe completely the behaviour of the 3-terminal component as 
measured at its terminals, we need to specify the values of only two of the 
current variables and two of the voltage variables. 

For a 2-terminal component, the graphical representation is a single edge 
with two variables associated with it (a current and a voltage). For a 
3-terminal component, we have four independent variables, so we need two 
edges for the graphical representation. One possibility is shown below: 


4a A 
2 
ae rN 
B Cc B c 
3-terminal component graphical representation 


We can express the variables i,, vp, etc. in terms of the edge variables v,, 
i, and v,, i, as follows, 


Clearly, the current flowing out of terminal B is equal to i,, so 
ig =i,. 
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fia 
LSS 
UAB UA 


| } 


Ls Je 
Bu ee Eee 


We show in Section 2 that Kirchhoff’s 
current law can be applied to any 
cutset in an oriented graph. 


We can assume that the component 
is connected to other components or 
to measuring instruments, so that 
there is such a cycle. 


In the graphical representation, the 
arrow on edge 1 points towards the 
vertex B, This means that the current 
i, is directed towards vertex B and 
therefore flows out of the 3-terminal 
component via terminal B — itis 
therefore equal to the current ig. 


Similarly, 

ic by 
and, since i, =—ig—ic, we have 

ig ts i, = h. 
The voltage across terminals A and B is clearly equal to v,, so 

Vag =. 
Similarly, 

U4 =~ Py 
and, since vgc = Vag - Uca, We have 

Vgc =U, - 2}. 
We showed earlier that a 2-terminal component is represented by one edge 
and has one component equation. We have just shown that a 
3-terminal component is represented by two edges and, by analogy with 
the 2-terminal component, has two component equations. As a result, we 
can regard any 3-terminal component as being composed of two 2-terminal 
components, one corresponding to each edge of the graphical 
representation. To represent the 3-terminal component shown 
diagrammatically above, we chose the two edges AB and AC. Instead, we 
could have chosen the two edges AB and BC, or the two edges BC and AC. 


In the following problem we ask you to find relationships between the 
variables for the different choices of pairs of edges. 


Problem 1.3 
The following graphs represent the same 3-terminal component: 
4 A 
yj \ 
a Dai ) a ae 
(a) (b) 


In each case, find i, i,, v, and v, in terms of the variables i,, i,, 0, and v, 
defined above. 


The method of representation just described can be extended to the 
graphical representation of n-terminal components for larger values of n. 
Just as a 3-terminal component can be represented by a tree with two edges, 
so an -terminal component can be represented by a tree with n — 1 edges. 
Moreover, we have seen that each of the labelled trees with three 
vertices (shown below) is a suitable representation of the 3-terminal 
component. 


B Cc B c B c 
Similarly, any labelled tree with n vertices is a suitable representation of 


an n-terminal component. Thus any n-terminal component can be regarded as 
n-1 2-terminal components, one corresponding to each edge of the tree. 


Problem 1.4 
(a) Find a graphical representation of the system shown in the margin. 
(b) How many component equations are associated with this system? 
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We conclude this subsection by listing some important types of component, 
together with their graphical representations and their corresponding 


component equations. 
Table 1.2 
graphical representation 4 _—- _8 
two-terminal through across component 
component variable variable equation 
resistance R 
$ 8 current i voltage v v=Ri 
resistor 
capacitance C 
4 current i voltage v i= 4? 
capacitor 
A__ inductance L r di 
ban current i voltage v v=La 
inductor 
Ae. 
é—_G—3 current i voltage v v=V 
voltage source 
<i 
$_a@p—3 current i voltage v i=] 
current source 
dampin 
coefficient c a 
B force f velocity u fecu 
damper 
mass m du 
é force f velocity u fomap 
mass 
A stiffness k ldf This relationship is often written in 
TD force f velocity u Une at the form f= kx, known as Hooke’s law, 
i where x is the extension of the spring 
spring due to an applied force f. 
valve ‘resistance’ R 
flow q pressure p p=Rq 


hydraulic valve 
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graphical representation \ £ 
B 
three-terminal through across component 
component variables variables equations 
mh 
currents voltages i 7 
i, and ip v, and v2 5 m, 
h= “iy it 
ideal transformer 
am nn teeth 
™\— ia 
to angular Th= hy Ty 
iT, he ond T, velocities s 
rn teeth o, @, and @, =~; 0 
2 


i 
displacements i= pf 
x and x2 L 
m=," 
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1.3. Dual electrical networks 2.4 


The idea of duality can be applied to (planar) electrical networks as 
follows. Given an electrical network, we draw a plane drawing of a 
corresponding oriented graph, form the dual of this planar graph as 
described in Graphs 3, and then construct the corresponding network, 
which we call the dual network. The method is illustrated below. 


A 4 

Cc aa 1 Bax Cc 
; : 

D 

given network oriented graph G 


x 
dual network dual oriented graph G* 


We start with a given network N and draw the oriented graph; the 
oriented graph G has four faces, which we have labelled w, x, y and z. We 
construct the dual oriented graph G* by drawing four vertices W, X, Y and 
Z corresponding to the four faces, and drawing an edge joining each pair of 
vertices separated by an edge in G. However, we are now dealing with 
oriented graphs, so we need to associate an arrow with each edge in G*. 
The rule for doing this is illustrated below. 


go clockwise from edge of G 
jo edge of G* in same direction 
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Finally, we construct the dual network N* from the dual oriented graph. 
It turns out that voltage in the original network N corresponds to current in 
the dual network N* and vice versa. Since these roles are interchanged, the 
roles of v and i in the component equations are interchanged, so we obtain 
the following correspondences between the variables and components in 
the two networks N and N*. 


dual network 
network N dual network N* 
voltage current 
current voltage 
resistor (resistance R) resistor (resistance 1/R) 
capacitor (capacitance C) inductor (inductance C) 
inductor (inductance L) capacitor (capacitance L) 


In Section 2 we find further correspondences between the fundamental 
equations for voltages and currents in the two networks. 


Problem 1.5 
Construct the dual of each of the following networks: 


electrical network oriented graph 
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14 Assumptions in modelling component behaviour 


In our discussion of component equations, we made a number of implicit 
assumptions about component behaviour. In particular, we assumed that 
the behaviour of a component does not change if it is moved from one 
system to another, and that the behaviour is time-invariant and linear. If 
these assumptions are not applicable to the components in hand, the 
mathematics will often turn out to be more complex and may require 
special manipulative techniques. We now comment briefly on the 
assumptions made in our treatment of physical networks. 


System independent behaviour 


One way of finding out how a component interacts with other components 
in a system is to remove it from the system and connect its terminals to a 
test system which can be used to take measurements. In this way we can 
construct a model of its behaviour, expressed as the component equation. 
There is, however, a very important implied assumption in doing this. We 
are assuming that the component will behave in the same way in our test 
system as in the original system. Otherwise, we would not be able to use 
the information obtained from the test to work out how the component 
behaves in the original system. Another way of putting this is to say that 
we must be careful that the model (or component equation) which we use is 
appropriate for the system in question. For example, one system may 
operate in a different type of environment from a second system, with a 
different ambient temperature, humidity or air pressure, and this may 
significantly affect the component behaviour. Even if the environments of 
two systems are the same, we may not be able to use the same model of a 
component for both systems. For example, a model of a transistor which is 
suitable for describing its behaviour in a hi-fi amplifier may not be 
suitable for describing its behaviour in a computer. 


Time-invariance 


The assumption that the behaviour of a component is time-invariant 
means that we assume its behaviour is the same today as it was yesterday 
and will be tomorrow. We must be aware, however, that the behaviour of 
real components will change with time. If this change is insignificant in 
the application we are concerned with, then the assumption of time- 
invariance is reasonable. 


Linearity 


Mathematically, the assumption of linearity means that the equations 
describing the behaviour of a component are linear equations. A system 
made up of linear components is itself linear. This means that we can use 
the many mathematical techniques that have been developed to handle 
linear equations to work out how the system responds to changes of its 
inputs. 

A linear system obeys the principle of superposition. This states that if the 
response of a system to an input x, is y,, and the response to an input x, is y,, 
then the response of the system when both inputs are applied together. is y, + y). 
Another way of looking at this is to consider the response of a system to a 
particular input, If the system is linear, and if, for example, the size of 
the input is doubled, then the resulting response will also be doubled. 


The following diagrams illustrate this principle for a simple electrical 
system in which the behaviour of the components can, for our purposes, be 
modelled by linear equations. In each case, a capacitor and a resistor are 
connected via a switch to a battery. 


1 volt 
t 2 volts 


battery 


capacitor 
(a) first case (b) second case 


In the first case, the voltage of the battery is 1 volt, and the system's 
response — the voltage v measured across the output terminals — on 
closing the switch is shown below, If we now replace the battery by one of 
2 volts and repeat the experiment, then the output voltage at any instant 
will be double what it was in the first case. 


voltage 


time 


2 Electrical networks: matrix equations 


In the previous section we showed that an electrical network consisting of 
a number of interconnected components can be represented by an oriented 
graph. We also discussed the equations relating the voltages and currents Remember that voltages are across 


associated with each component — the component equations. variables and currents are through 
iables. 
In this section we discuss two further sets of equations — the equations eee 


obtained by applying Kirchhoff’s voltage law, and the equations 


obtained by applying Kirchhoff’s current law to the oriented graph. In 
Section 3 we show that these three sets of equations: 


* the component equations, 
*  Kirchhoff’s voltage law equations, 
*  Kirchhoff’s current law equations, 


are sufficient to enable us to solve electrical network problems — that is, 
to find the values of the through and across variables associated with all 
the components of a given network. 


We can apply Kirchhoff’s voltage law to every cycle of the oriented 
graph, and we can apply Kirchhoff’s current law to every vertex of the 
graph. However, the equations which result contain a good deal of 
redundant information. What we should like is a systematic method for 
obtaining just enough equations to enable us to solve electrical network 
problems — or, to put it another way, to obtain sets of equations which 
contain all the information which Kirchhoff’s laws can tell us about the 
network, without any redundancy. 


Such a systematic method was given by Kirchhoff in 1847 in a paper 

entitled ‘On the solution of the equations obtained from the investigation 

of the linear distribution of galvanic currents’. His method used the We introduced the concept of a 
concept of what is now called a spanning tree and involved finding a set of _ spanning tree in Graphs 2. 
fundamental cycles and a set of fundamental cutsets. We use the same 

method, although the language in which we express it is somewhat 

different from that used by Kirchhoff. We consider only connected graphs 

— this means that we can always find a spanning tree. 


Our aim is to describe a method of solving network problems which is 
appropriate for computer analysis, so we need to overcome the problem of 
obtaining the fundamental cycles and cutsets using a computer. We 
describe a method which enables us to obtain the fundamental cycles and 
cutsets directly from the incidence matrix of the oriented graph. The great 
advantage of this method is that the only information about the oriented 
graph which we have to feed into the computer is the incidence matrix 
itself, and we can obtain the incidence matrix by inspection from the 
graph. 


2.1 Kirchhoff’s voltage law equations 


We can apply Kirchhoff’s voltage law to any cycle of the oriented graph 
of an electrical network. But if we apply this law to all the cycles, we end 
up with a set of equations which contain a lot of redundant information. 
We demonstrate this for the following electrical network which we use as 
an illustrative example throughout this section. 


i 3 ¥ Aa? 6 § 
oe xt it te s A | 
F 4 E joes D F ri. E oF D 
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electrical network oriented graph 
Let us start by applying the voltage law to the cycles BCEB, CDEC and _ Recall that in going around each cycle 
BCDEB of the oriented graph. we take a voltage to be positive if it is 
hi ; , oriented in the same direction as the 

This gives the following voltage equations. cycle, and negative if it is not. 

BCEB 5 +07 —05=0 

CDEC 05 +0 — 0 =0 

BCDEB U6 +09 +0 —¥5 =0 
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The last of these equations is simply the sum of the other two, and so gives 
no further information about the voltages. 


Similarly, applying the voltage law to the cycles ABEFA, BCEB, CDEC 
and ABCDEFA (using the right-hand edge AF), we obtain the following 
voltage equations. 


ABEFA +05 +%4—2%=0 
BCEB U5 +07 —V5=0 
CDEC Wy +0,—07=0 
ABCDEFA D3 +05 + Vy + Uy +14—%=0 


Again, the last of these equations is the sum of the others, and so gives no 
further information about the voltages. 


In each of these cases, we can write each equation as the sum or difference 
of the other equations, and we express this by saying that the equations 
are linearly dependent. More generally, a number of equations are said to be 
linearly dependent if at least one of them ‘depends’ on the others, in the 
sense that it can be obtained by adding or subtracting multiples of the 
others. If none of the equations depends on the others in this way, then the 
equations are said to be linearly independent. For example, the voltage 
equations for the cycles ABEFA, BCEB and CDEC are linearly 
independent, since each contains a term contained in neither of the other 
two equations. 


Problem 2.1 


Consider the following oriented graph, which corresponds to the network 
in Problem 1,1. 


electrical network oriented graph 


(a) Write down the voltage equation for each of the seven cycles of the 
graph. 


(b) What is the maximum number of these seven equations which are 
linearly independent? 


2.2 Fundamental cycles 


We now investigate the problem of finding a set containing the maximum 
number of Kirchhoff’s voltage law equations which are linearly 
independent. What we need to do is to find a set of cycles in an oriented 
graph which give rise to such a set of linearly independent equations. We 
say that a number of cycles in a graph are linearly independent if the 
corresponding voltage equations are linearly independent. For instance, in 
our illustrative example, the cycles ABEFA, BCEB and CDEC are linearly 
independent, since they correspond to linearly independent equations. It 
follows that, in order to find a set of linearly independent equations for 
the network, it is enough to identify the corresponding linearly 
independent cycles. We can find such a set of cycles quite easily; to do so, 
we use Kirchhoff’s idea of a spanning tree. 


Definitions 


Let G be a connected graph. A spanning tree in G is a subgraph of G 
which includes all the vertices of G and is also a tree. The edges of the 
tree are called branches. If we remove from G all the edges of a 
spanning tree, the remaining subgraph is called a co-tree and its edges 
are called chords. 


The following diagram shows the graph in our example, one of its 
spanning trees, and the associated co-tree, 


A B c A B c 
branches chords 
F E D F EB D F E D 
graph spanning tree co-tree 


For any connected graph G, we can find a spanning tree systematically. 
We start by choosing any cycle in G and removing one of its edges. 


Since we cannot disconnect a graph by removing just one edge from a cycle, 
we still have a connected graph. We now repeat the process until there 
are no cycles left, and there is our spanning tree. For the above graph, 
removal of the edge CD from the cycle CDEC, the edge BC from the cycle 
BCEB, the edge BE from the cycle ABEFA, and one of the edges FA from 
the cycle AFA, leaves the spanning tree shown. 


4 3 


Problem 2.2 


(a) Find two more spanning trees in the above graph, and list the 
branches and chords in each case. 


(b) Ifa connected graph has n vertices and m edges, how many branches 
and chords are there in each spanning tree? ) 


(c) What happens if we add a chord to a spanning tree? Illustrate your 
answer using the spanning tree given in our example. 


The last part of Problem 2.2 shows us how to find a set of linearly 
independent cycles. We first choose a spanning tree; then we add one of 
the chords to the tree — this produces exactly one cycle. We repeat this 
procedure for each chord. If the graph has k chords of the spanning tree, 
then this procedure produces k cycles. Since each cycle contains one chord 
which is not contained in any of the other cycles, the resulting set of cycles 
must be linearly independent. 


If we were seeking a set of equations to solve by hand, we would choose the 
spanning tree with some care. In particular, we would choose the tree so 
that in general each edge for which the across variable was known would 
be a branch of the tree, and each edge for which the through variable was 
known would be a chord. However, this selection is not necessary in the 
mathematical sense; it serves only to make the resulting equations quicker 
to solve by hand. We shall describe a method appropriate for solving 
physical network problems on a computer, where such considerations are 
unimportant, and so we do not need to follow a selection procedure: we 
simply choose any convenient spanning tree. 
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‘The number of spanning trees in a 
graph can be very large; for instance, 
the graph in our example (shown 
below) has exactly 50 different 
Spanning trees. 


Note that a co-tree is, in general, not a 
tree. 


If there are no cycles, the graph Gis 
itself a spanning tree. 


Since the set of cycles we have found using our spanning tree are linearly 
independent, the corresponding voltage law equations must also be 
linearly independent. This leads to the following definitions. 


Definitions 


Let G be a connected graph, and let T be a spanning tree of G. The set of 
fundamental cycles associated with T consists of the cycles of G each of 
which is obtained by adding a single chord to T. The corresponding 
voltage law equations are called the fundamental cycle equations. 


For example, for the spanning tree shown below 


fes 4 4 B 
branches chords. 
F E D FP E D F E D 
graph spanning tree co-tree 
we have the following table: 
fundamental fundamental 
chord cycle diagram cycle equation 


BC ABCEFA D3 +05 +d, + 04-1)=0 
BE ABEFA ty +05+ 4-1 =0 
cD CDEC V4 9+ Ug 0, =0 
E 8 D 
A 
FA AFA’ 1 +0, =0 
F 


It follows from Problem 2.2(b) that, if G has n vertices and m edges, then 
any spanning tree gives rise to m —n +1 chords, and so there are m—n + 1 
cycles in a fundamental set. This number is the largest number of linearly 
independent cycles in G, and is called the cycle rank of G. It follows that 
there are m~—n + 1 fundamental cycle equations, and that this is the 
largest possible number of linearly independent cycle equations. We 
conclude that, in any electrical network problem, we can obtain all the 
information contained in Kirchhoff’s voltage law by: 


(a) choosing a spanning tree; 
(b) finding the fundamental cycles associated with it; 
(c) finding the corresponding fundamental cycle equations. 


Note that the cycles are taken in the 
same direction as the chords they 
contain, Thus we have written ABEFA 
and CDEC rather than AFEBA and 
CEDC, G 


An analogous result holds for the 
information contained in the 
generalized form of Kirchhoff’s 
voltage law for other types of network. 
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Problem 2.3 
For each of the spanning trees shown below, write down: 
* a list of the chords; 

* the corresponding fundamental cycles; 

* the corresponding fundamental cycle equations. 


In each case, show how to obtain the voltage law equation for the cycle 
ABCDEA from the fundamental cycle equations. 


Before leaving fundamental cycles, we show how the fundamental cycle 
equations can be represented in matrix form. The rows of the matrix 
correspond to the chords associated with a spanning tree, and the columns 
correspond to all the edges of the graph G. If G has n vertices and m 
edges, we get an (m—n +1) xm matrix. This matrix is called the 
fundamental cycle matrix, and is denoted by Cy. 


To construct this matrix, we add each chord to the spanning tree in turn, 
and look at the edges in the corresponding fundamental cycle. Tracing 
around this cycle in the direction of the chord, we fill in the row 
corresponding to the given chord as follows: we write 


1 in each column corresponding to an edge of the cycle oriented in the 
same direction as the chord; 


-1  ineach column corresponding to an edge oriented in the opposite 
direction; 


0 in each column corresponding to an edge which is not in the cycle. 
For example, for the fundamental cycle ABCEFA illustrated below, the 


row corresponding to the chord BC has 1 in the columns corresponding to 
AB, CE, EF and BC, -1 in the column corresponding to AF, and 0 elsewhere. 


Repeating this procedure for all four fundamental cycles, we obtain the 
following matrix: 


AB AF CE DE EF BC BE CD FA 


Bef acl 1) <0) “dy aio nomen 

BE| d=1, 0.10). Oaeae Deo 

€D\ 0) (0 =k 1 0 derolete Dilip 

FA 0 1 0) 0 0 On Omori 
branches chords 


Having found the fundamental cycle matrix Cs, we can write the 
fundamental cycle equations in matrix form as 


Crv=0, 
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It is usual to write the columns. 
corresponding to the branches of the 
spanning tree first, followed by those 
for the chords. 


where v is the column vector of edge voltages, written in the same order as 
the columns of the matrix, and 0 is the appropriate zero vector. For 
example, the above matrix gives 


edgenumber > 3 2 7 8 4 6 5 9 1 
Pt 0 4 60) 102 01s fo 
1 <1 0 0: {0 1 0 O}leI_h0 
60 0-2 10 0 O 1 Oj/y| Io 
0 10 0 0 0 O O Tila} 10 
U% 
% 
Us 
U9 
y 


Multiplying the two matrices on the left, we obtain the following 
equation, which is equivalent to the four equations on page 27. 


03 — U2 +07 +04 +06 0 
03-02 +04 +05 0 
07 +g +09 0 

U2 +, 0 


Problem 2.4 


For each of the spanning trees in Problem 2.3, find the corresponding 
fundamental cycle matrix, and use it to derive the fundamental cycle 
equations in each case. 


‘2.3 Kirchhoff’s current law equations 


We now turn our attention to Kirchhoff’s current law, which states that 
the algebraic sum of the currents at each vertex of the network is zero. As 
with the voltage law, we find that if we apply the current law to every 
vertex of a graph, then we end up with more equations than we need. 


For example, consider our oriented graph: 


Applying the current law at each vertex, we obtain the following 
equations. 


vertex A }-i-15=0 
vertex B is —i 

vertex C 

vertex D 

vertex E ist+i7+ig—ig=0 
vertex F int+ig-i, =0 


These six equations are not linearly independent: if we add any five of We show later that for a graph with n 
them, we get the sixth. However, any five of these equations are linearly _ vertices there are always exactly n-1 


independent. 


linearly independent equations. 
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In fact, for the solution procedures which we use here, it is more 
appropriate to apply Kirchhoff’s current law to the cutsets of the graph, 
rather than to the vertices. 


To see what we mean by this, consider the cutset consisting of the edges 
BC, CE and DE of our graph: 


Removal of the edges of this cutset disconnects the graph into two parts, 
X and Y. If we add the equations for the vertices in just one of these parts 
(Y, say), we find that the currents along those edges lying entirely within 
that part cancel out, and we are left with an equation relating the currents 
in the edges of the cutset, 


For example, if we add together the above equations for the vertices C and 
D, the current ig cancels out, leaving the equation 


Ebb = 0, 
This equation relates the currents in the cutset edges BC, CE and DE. Note 


that we get the same equation if we add the equations for the vertices A, 
B, Eand F in the set X. 


We can generalize this idea, and use the current law to find an equation 
relating the currents in any cutset. If the cutset separates the graph into 
two parts X and Y, then the equation has the form 


(the sum of the currents from X to Y) 
—(the sum of the currents from Y to X) 
=0. 
Or simply, 
the algebraic sum of the currents in any cutset is zero, 


However, if we repeat this procedure for all the cutsets in the graph, we 
find that the resulting equations are not all linearly independent. 


2.4 Fundamental cutsets 


As with the cycle equations, we wish to find the largest possible number of 
linearly independent cutset equations. To do this, we find a set of 
fundamental cutsets. As with fundamental cycles, these are constructed 
using a spanning tree. 

We illustrate the method with the graph and spanning tree we used 
previously: 


The removal of any branch separates the tree into two parts, X and Y; for 
example, the removal of the branch EF separates the tree into two parts, 
one containing the vertices C, D and E, and the other containing the 
vertices A, B and F. Since the branch EF is oriented from E to F, we let X be 
the part containing E, and Y be the part containing F. We now list all the 
edges of the oriented graph joining a vertex in X and a vertex in Y — these 
are the edges EF, BE and BC. One of these edges is the branch EF that we 
removed to disconnect the tree; the others are the chords BE and BC; 
together they form a cutset of the network graph. 


We repeat this procedure for each branch of the spanning tree. The 
resulting cutsets are called fundamental cutsets; since there are n - 1 
branches in the spanning tree, there are n — 1 fundamental cutsets. This 
leads to the following definitions. 
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We defined a cutset of a connected 
graph in Networks 1. 


The word ‘algebraic’ means that we 
take account of the orientation of 
each edge, 


Definitions 


Let G be a connected graph, and let T be a spanning tree of G. The set of 

fundamental cutsets associated with T consists of the cutsets of G 

obtained by removing a branch of T, thus separating the tree into two 

parts, X and Y, and listing the edges of G joining a vertex in X and a 

vertex in Y. The corresponding current law equations are called the 
fundamental cutset equations. 


We illustrate this definition by our example, using the same spanning tree 
as before. 


vertices vertices fundamental 
branch = inX inY fundamental cutset diagram cutset equation 
AB A’C,DE,F B {AB, BC, BE} i, ~is-i, =0 
AF A, B C,:D, E, F {AF, BC, BE, FA} itis +i,-i,=0 
CE c A,B, D,E,F (BC, CD, CE} Zl ip +in-ig=0 

F E D 

DE D A,B,C, E,F (CD, DE} ig-iy =0 
EF CG, DE A, B, F {BC, BE, EF} 


F E 


a 
' 
Ya 


x 


We say that a number of cutsets in a graph G are linearly independent if 
the corresponding current law equations are linearly independent. Thus 
the cutsets in a set of fundamental cutsets are all linearly independent, 
since each cutset contains one edge (the branch of the spanning tree) not 
contained in any of the other cutsets. If G has n vertices, then any 
spanning tree has n — 1 branches, and so there are n — 1 cutsets in a 
fundamental set. This number is the largest number of independent cutsets 
in G, and is called the cutset rank of G. It follows that there are n -1 
fundamental cutset equations, and that this is the largest possible number 
of linearly independent cutset equations. Moreover, these equations, 
rather than the equations obtained by applying Kirchhoff’s current law 
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at the vertices of G, are the appropriate equations for our purposes. We 
conclude that, in any electrical network problem, we can obtain all the An analogous result holds for the 
information contained in Kirchhoff’s current law by: information contained in the 

i i : generalized form of Kirchhoff’s 
* (choosing a.spanning tree; current law for other types of network. 
* finding the fundamental cutsets associated with it; 


* finding the corresponding fundamental cutset equations. 


Problem 2.5 
For each of the spanning trees shown below, write down: 
(a) a list of the branches; 

(b) the corresponding sets X and Y; 

(c) the corresponding fundamental cutsets; 

(d) the corresponding fundamental cutset equations. 


In each case, show how to obtain the current law equation for the cutset 
{AB, BC, BD} from the fundamental cutset equations. 


As with fundamental cycles, we can represent the fundamental cutset 

equations in matrix form. The rows of this matrix correspond to the 

branches of a spanning tree, and the columns correspond to all the edges of As we said earlier, it is usual to write 
the graph G. If G has n vertices and m edges, we get an (n —1)x m matrix. the columns corresponding to the 


This matrix is called the fundamental cutset matrix, and is denoted branches first, followed by those for 
by Dy. the chords. 


To construct this matrix, we take each branch of the spanning tree in turn, 
and look at the edges in the corresponding fundamental cutset. We then 
fill in the row corresponding to the given branch as follows: we write 


1 in each column corresponding to an edge of the cutset oriented 
from X to Y; 


-1 in each column corresponding to an edge of the cutset oriented 
from Y to X; 


0 ineach column corresponding to an edge not in the cutset. 


For example, for the fundamental cutset (AB, BC, BE} illustrated in the 
margin, the row corresponding to the branch AB has 1 in the column 
corresponding to AB, ~1 in the columns corresponding to BC and BE, and 0 
elsewhere. 

Repeating this procedure for all five fundamental cutsets (shown on 
page 31), we obtain the following matrix: 


AB AF CE DE EF BC BE CD FA 


aBaf 10 0 0 0-1-10 0 
af 0! % OO) by BO 10) 4, 
CED) 0 lo 10) SO aos 0: branches 
DE} O 0 0 10 00-1 0 
BPO) 100 i 07 to T onl Oa 
branches chords 
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Having found the fundamental cutset matrix D;, we can write the 

fundamental cutset equations in matrix form as 
Dyi=0, 

where i is the column vector of edge currents, written in the same order as 

the columns of the matrix. For example, the above matrix gives 


edge number ~ 312 2 68 Cee tae: hee: Jem 
10 0 0 0-1-1 0 O}f%] fo 
Oi 2-0 (0) 0) EO) ane Te: 
00 20 0-4 0 1 Ole, Jel0 
0 0 0 1 0 0 0 -1 O}fig} JO 
0060 0 1-1-1 0 OG) [0 


Multiplying the two matrices on the left, we obtain the following 
equation, which is equivalent to the five equations on page 31. 
is-is-is] [0 
i, +igtis—i,] |0 
ip ig tig |=|0 
0 
0 


9 
bias 


Problem 2.6 


For each of the spanning trees in Problem 2.5, find the corresponding 
fundamental cutset matrix, and use it to derive the fundamental cutset 
equations. 


2.5 Obtaining the fundamental cycle and cutset matrices 


So far we have shown that: 

* all the information which can be obtained by applying Kirchhoff’s 
voltage law to an electrical network is contained in the matrix 
equation 

Crv =0, 
where C; is the fundamental cycle matrix; 


* all the information which can be obtained by applying Kirchhoff’s 
current law to an electrical network is contained in the matrix 
equation 

Dri=0, 
where Dy is the fundamental cutset matrix. 


We have shown how the fundamental cycle matrix C; and fundamental 
cutset matrix Ds with respect to a given spanning tree can be obtained by 
inspection of the oriented graph of an electrical network, Even for the 
small networks we have considered, this is a somewhat tedious process, so 
we would like to be able to get a computer to do this for us. The problem 
with this is that the method we have described is not very suitable for 
adaptation as a computer program. Fortunately, however, the matrices 
C; and Dy can be obtained directly from the incidence matrix of the 
oriented graph by an elegant method which is suitable for translation into 
a computer program. We now describe this method, beginning with the 
definition of the incidence matrix of an oriented graph. 


The incidence matrix of an oriented graph 


We define the incidence matrix of an oriented graph by regarding the 
oriented graph as a digraph in which the direction of each arc is just the 
reference direction of the corresponding edge of the oriented graph. 
Adapting the definition of the incidence matrix of a digraph given in 
Graphs 1, we obtain the following definition. 


Definition 


Let G be an oriented graph with n vertices and m edges. The incidence 
matrix B(G) is the n x m matrix in which the entry in row i and column j 
is 


1 if edge j is incident with, and oriented away from, vertex i; 
-1_ if edge; is incident with, and oriented towards, vertex i; 


0 if edge j is not incident with vertex i. 


For our example, the incidence matrix is 

AB AF BC BE CD CE DE EF FA 
alt ft 0) O° 10) 10) 60) 0 =1 
Blete 10 10s Oe 100 0 x0: 
tO 30) ee 0h sa ON 0 60! 
pho Oo O Q1 0 1 OO 
E}O- 0 O-1 O-1 -1 1 0 
FIO) =1 (0: 10) 90) 0) Ol = 2 


Each column of the incidence matrix contains one 1 and one -1, so that if 
the six rows are added together, the result is a row of zeros. Thus each row 
is equal to minus the sum of all the other rows. This result holds for the 
incidence matrix of any oriented graph. Since the rows of a matrix are 
linearly dependent if one of them can be obtained by adding or subtracting 
multiples of the others, the rows of any incidence matrix are linearly 
dependent. It follows that we can omit any single row of B(G) without loss 
of information. For example, in the above matrix, if we omit the row 
corresponding to the vertex C, then we get the following matrix. We 
denote it by Bo: 


AB AF BC BE CD CE DE EF FA 


Ay kA) 50) AO OG) OeiOi ed: 
Be OPER MeO GO 102s 0) 
Dy Ol O° 300 0) =A 50) 3) 10) 30.) Bo, 
eO 0 0-1 0-1-1 1 °0 
F100) —E SOON OIG AO 


Such a matrix is called a reduced incidence matrix, and the vertex C is 
called the reference vertex. 


If we are given the matrix Bg, we can recover the original incidence matrix 
B(G) by inserting the missing row (calculated by taking minus the sum of 
the rows of By). However, we cannot replace two missing rows, since we do 
not know which direction to assign to an edge joining the corresponding two 
vertices. It follows that if G has n vertices, then the maximum number of 
linearly independent rows of the incidence matrix B(G) is n — 1. 


Suppose that we have chosen the spanning tree of G with branches AB, 
AF, CE, DE and EF. We can rearrange the columns of the reduced incidence 
matrix By so that those columns corresponding to branches of the tree come 
first, followed by those columns corresponding to the chords, as follows. 
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apt no: 0) OO: 0 0 =r 
aj-1. 0 0 0 O11 1 0 O 
Bo=D/0 0 0 1 010 0-1 O 
FO 0-1-1 110-12 0 Oo 
FLO -1 0 0-110 0 0 1 
branches chords 


We have partitioned Bo into two parts, one part corresponding to the 
branches, and one part corresponding to the chords. We can abbreviate 
this partitioned matrix by writing Bp = [B, | B.], where B, (the tree part) 
consists of the first five columns, and B, (the co-tree part) consists of the 
last four columns: 


AB AF CE DE EF BC BE CD FA 
At i 10) 0) 10 a4fo 0 0 -1 
Bi-l 0 0 0 O Bi rt 0 0 
B=2}0 0 0 1 O|andB,=D) 0 0 - 0 
e}O 0-1-1 1 BO -1 0 0 
FLO -1 0 O -1 FLO 0 0 1 
Problem 2.7 


Consider the oriented graph and spanning tree shown in the margin. 
Write down: 

(a) the incidence matrix B of G; 

(b) the reduced incidence matrix By, using vertex D as reference vertex; 
(c) the matrices B, and B,, using the spanning tree shown. 


Obtaining the fundamental cycle and cutset 
matrices from the incidence matrix 


The reduced incidence matrix By completely specifies the corresponding 
oriented graph — that is, it tells us which edges join which vertices, and 
it tells us the reference direction of each edge. It is a convenient 
representation of the graph to use for computer analysis, since it can easily 
be written down by inspection of the oriented graph, and it is suitable for 
feeding into a computer. Because the reduced incidence matrix completely 
describes the oriented graph, it should theoretically be possible, once a 
spanning tree has been specified, for the computer to generate the 
fundamental cutset and cycle matrices without further reference to the 
graph itself. We now describe how this can be done. 


First we must specify a spanning tree. Any spanning tree will serve our 
purpose, and an appropriate algorithm (such as the one given on page 24) 
can be used to find a spanning tree using the information contained in the 
reduced incidence matrix. We can then arrange the reduced incidence 
matrix in the form By = [By | B,], as described above. 


Suppose that the oriented graph G has n vertices; then By has n -1 
linearly independent rows (corresponding to the n — 1 linearly independent 
rows of the reduced incidence matrix) and n — 1 columns (corresponding to 
the n — 1 edges in the spanning tree). So B, is a square matrix whose rows 
are linearly independent. It follows from the theory of matrices that By 
has an inverse matrix; that is, there exists a matrix By! of order (n — 1) x 
(n-1) such that 
Be 1B, = By Bet = Ty -1 
where I, _; is the identity matrix of order (n — 1) x (n - 1). 
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We now state (without 
matrices Cr and D; from 


roof) a result which enables us to obtain the 
e incidence matrix. 


Theorem 2.1 


The fundamental cycle matrix Cy and the fundamental cutset matrix 
D; can be expressed in terms of B; and B, as follows: 


C= [Be BI —n +1) 
Dp= [In -1! Be *Be]. 


fundamental cycle matrix 


fundamental cutset matrix 


We shall explain in detail how these expressions are used to calculate 
C; and Dy, illustrating the procedure using our example. Unfortunately, 
the proofs of the validity of these expressions for Cr and Dy are too 
complicated to be included here. 


We first turn to the problem of finding By". In general, it is not easy to find 
an inverse matrix. However, for this particular type of matrix B,, there is 
a simple method for finding its inverse matrix By !. 


STEP 1 We first draw the spanning tree, marking in the orientations of 
the edges, and indicating the reference vertex. 


STEP 2 We next trace the path from each vertex of the tree to the 
reference vertex, noting carefully which branches are traversed in the 
same direction as the reference direction, and which branches are 
traversed in the opposite direction — for example, the path from vertex A 
to the reference vertex C consists of the branches 


AF (forward), EF (backward), CE (backward), 
STEP 3 Finally, we construct the matrix whose columns correspond to the 
vertices (other than the reference vertex) and whose rows correspond to 
the branches of the spanning tree. Each column represents the path from 


the corresponding vertex to the reference vertex; the element appearing in 
the column is 


1 _ if the corresponding branch is traversed in the same direction as 
the reference direction; 


-1__ if the branch is traversed in the opposite direction; 
0 if the branch is not included in the path. 


For the above example, we obtain the following matrix. 


A 8 DS F 


an) O -1 0 0 O 
agi 2 61. 10 60 0 
cE}-1 -1 -1 -1 -1 
be OO 0 61 UO C0 
EFi/-1 -1 0 O -1 


This is the required matrix By, We can check this by forming the matrix 
products B,'B, and B,B,7: 


O=1 0) 0 Of 1 io oMto7N (noo ome 
2 1 1D) 0 (Oi S0 10 rae SOON NOt cao) 0 

BB, =}-1 -1 -1 -1 -1]]0 0 0 1 O|=/0 0 10 O|=1; 
0 0 1 0 Of/o 0-1-1 1] jooo10 
-1-1 0 0-1//0 1 0 0-1} |joo001 
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1p.) is the transpose of By 'Be. 
t pos 


4 B reference 
= vertex 
| gO 
F E B 
£ 


TL es CON Ve iar 1" 50)) 10" 10: 10000 
S20) 0) HOO) Lod ONO) 0/0 Ct s0:50) 0 
BB'=| 0 0 O 1 O|/-1 - -1 -1 -1/=]0 0 10 O|=1 
Ot UO! Gil © Die} 0) 00nd 
O -1 0 0 -1)/-1 -1 O 0 =| jo O00 01 
Problem 2.8 


Use the above method to find the inverse B,? of the matrix B, in Problem 
27, and check your answer by calculating By'B, and B,By". 


Having constructed B,', we can use it to calculate the fundamental cutset 
and cycle matrices using the expressions given in Theorem 2.1. We 
illustrate this procedure by applying it to our example. 


First we form the matrix product B,'B,: 


BC BE CD FA 

alO -1 0 0 O]f0 0 0 -17 

ABT 10) GO Willd 1) cond 

B,"B, =CE/-1 -1 -1 -1 -1/]0 0 -1 0 
pb} OO 60 61 (0 «CO}} 0 -1 OO CO 
FFl-1-1 0 0 -1/[/0 0 0 1 


BC BE CD FA 
ABl=1 = 0) 10 
AF} 1 1 -0 =1 

=cel-1 0 1 0 
De} Oo 0 -1 0 
EFi-1 -1 0 0 


We now insert the rows and columns of the identity matrix I; at the left- 
hand side of the matrix: 


AB AF CE DE EF BC BE CD FA 

afl 0 0 0 O!-1 -1 0 0] 

a}0 1 0 0 011 1 0-1 
[is|BB.J=ce}0 0 1 0 0/1 0 1 0 
p00 0 1 010 0-1 0 

FLO 0 0 O 1;-1 -1 0 Of 


This is the fundamental cutset matrix Ds, as can be seen by comparing it 
with the matrix obtained on page 32 by inspection of the oriented graph. 


To form the fundamental cycle matrix C;, we take the transpose of By"B, : 


AB AF CE DE EF 

peal 1 <1) 0 1 
Bectepe|ed dt 10) 0" a 
(8B) =)0 0 140 
ral 0-1 0 0 0 


We then multiply each element of this matrix by -1, and write the 
identity matrix I, to the right of it: 


AB AF CE DE EF BC BE CD FA 

Laps) SUB Wn sel Gh Gh 
Savile tat 0 0 110 1:0 0 
[8 Be) In]=25 00-1 10!0010 
flO 10 00/0001 


This is just the fundamental cycle matrix C,, as can be seen by comparing it 
with the matrix we obtained on page 28 by inspection of the oriented 
graph, 
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Problem 2.9 


(a) Use the matrices B, and B,? calculated in Problems 2.7 and 2.8 to 
find the corresponding fundamental cutset and cycle matrices Dy and 
C;, and check that your answers agree with those of Problems 2.6 and 
2.4 respectively. 


(b) Calculate the matrix product Cpe. 


Problem 2.10 


For each of the following oriented graphs, find the fundamental cycle 
equations and the fundamental cutset equations, and comment on your 
results. In each case, use the spanning tree indicated by thick edges. 


(a) (b) 


2.6 Tellegen’s theorem ‘XU. 


We conclude this section by showing how to prove a result of great 
theoretical importance, known as Tellegen’s theorem, by making use of the 
incidence matrix. It may be stated in the following form. 


Theorem 2.2: Tellegen’s theorem 


Suppose that there are two electrical networks which can be 
represented by the same oriented graph with n vertices and m edges. 
Then 


m m 
Lorik = 0, Logit = ply + dai +... + Onin 
kel 1 

where vx is the voltage associated with the kth edge for the first 


network, and i, is the current associated with the kth edge for the 
second network. 


Proof 


We start by applying Kirchhoff’s current law to each of the n vertices of 
the oriented graph. It is not difficult to see that the resulting equations 
can be expressed in the matrix form: 


BJ =0, 
where 
Bis the n x m incidence matrix of the oriented graph; 
J is the mx 1 column vector of currents in the second network 
associated with the m edges of the graph; 
0 is the nx 1 zero vector. 
We now turn to Kirchhoff’s voltage law. Since the voltages associated 


with the edges of the graph satisfy Kirchhoff’s law, we can assign a 
unique potential to each of the n vertices such that the voltage v, in the 


first network associated with edge k is given by 
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1 = P,-Py 
where P, and P, are the potentials at the two vertices joined by the kth 


edge. The equation analogous to the above for all the vertices of the graph 
can be expressed in the matrix form: 


B'P=V, 
where 
BT is the m x n transpose of the incidence matrix B; 
P__ is the n x 1 column vector of potentials at the n vertices; 
V__ isthe mx 1 column vector of voltages (or potential differences) 
associated with the m edges of the graph. 


We can write the expression in Tellegen’s theorem in matrix form thus: 


that the expression on the left is a 


mn a The brackets are inserted to indicate 
Lreik |= VT, 
k=l 1x 1 matrix, 


where V7 has order 1 x m, and J has order m x 1, so V"J has order 1 x 1. 


We know that V = B! P, so we substitute this in the above matrix 
expression and thus obtain: 


m 
[ex] =(B" Py} 
=(P7B)J (BT P)" =P"(B")"=P"B) 
=P(BJ) (property of matrix multiplication) 
=PQ. x1 (BJ =O, x1) 
=001 (PT has order 1 x n). 


m 
Hence Sri, = 0, as required. a 
fet 


Tellegen’s theorem holds for any two networks which have the same 
oriented graph. Equally, we can apply it to just one network, in which case 
the result is that Z ogi, = 0, where vj and ij are the voltage and current 


es 

associated with the Hin edge of this network. The product v,ix is the 
power (or rate of change of energy) associated with the kth edge, so the 
equation £ vi = 0 tells us that energy is conserved in the network. This is the law of conservation of 

all edges energy for a network. 
Note that we have derived this well-known law of physics for an Although we have stated Tellegen’s 
electrical network from Kirchhoff's current and voltage laws. In fact, theorem for electrical networks, it also 
given any two of these three laws (the conservation of energy law, holds for the through and across 


Kirchhoff’s current law and Kirchhoff’s voltage law), it is possible to Variables of any physical systems 
derive the third. 8 ) P which obey generalized forms of 
. Kirchhoff’s laws. 


Summary 


On the next two pages, we summarize what we have learned about the 
two examples of this section. These results will be used in the next section. 
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Example in text 


network 


oriented graph 


4 B iS 
spanning tree 
F E D 
A B 4 B A 
fundamental cycles wa )] ‘ 
F E F E E D F 
fundamental cycle equations Wty +U,+%-v,=0, 0, +5+%,-v,=0, U+t%-v,=0, v,+v,=0 
AB AF CE DE EF BC BE CD FA 
Bef Dot 80) i as Oe Og: 
‘ BE} 1-1 0 Of 1140 1 0 0 
fundamental cycle matrix Crecolo oO <1 1 0 NiO wadinet chords 
Alo 1 6 0 of 0 0 O i 
branches chords 


fundamental cutsets 


fundamental cutset equations 
AB AF CE DE EF BC BE CD FA 
48,31 0 0 O Oj-1 -1 0 O 
AF} O 1 0 0 0; 1 10 -1 
fundamental cutset matrix D,=cE} 0 0 1 0 O1-1 0 1 0}> branches 
DE} O 0 0 1 OO O-1 O 
FLO 0 0 O i}-1 -1 0 O 
branches chords 


Example in problems 


network 


oriented graph 


4 
spanning tree : te 
D 
E fa 4 B 
fundamental cycles Cc a cel 
D E ic 
D 
fundamental cycle equations Uy -%%-0;=0, 0, +0,-v5,=0, v, +0, +0; +0, =0 
AB BC CD CE BDDEEA 
: Bo} 0 -1 -1 0;100 
fundamental cycle matrix Cr=DeE| 0 O 1 -110 1 O}} chords 
EA Ue) Ou I 
Se 
branches chords 
fundamental cutsets {AB, EA}, {BC, BD, EA}, (BD, CD, DE}, (CE, DE, EA} 
fundamental cutset equations i,-i, =0, i, +i,-i,=0, ig+i,-i, =0, 
AB BCCD CE BD DE EA 
ABiA (0! 0" 0)| 10) “0! =I 
fundamental cutset matrix Bc]0 10 0} nail 
eplo 0 10} 1 «1 off branches 
celO 00 1!0 1-1 


branches chords 
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‘3 Electrical networks: 
solving the network equations 


In the previous two sections we obtained a number of equations relating the 

voltages and currents in an electrical network. These equations are of three 

types: 

1. component equations, relating the voltages and currents in each 
component; 


2. fundamental cycle equations, relating the voltages in each 
fundamental cycle of the network; 


3. fundamental cutset equations, relating the currents in each 
fundamental cutset of the network, 


These three sets of equations contain all the information we need to be 
able to work out all the voltages and currents in the network. If the 
oriented graph of the network has n vertices and m edges, then there are 
m voltages and m currents — a total of 2m unknowns related by 


1. mcomponent equations; 
2. m-—n-+1 fundamental cycle equations; 
3. n~1 fundamental cutset equations. 


There are, therefore, m + (m—n +1) + (n -1) = 2m equations in 2m 


unknowns, 


In this section we show how these 2m equations can be organized 
systematically into a standard matrix form. We look first at networks 
that contain only resistors and independent sources, and show that, given 
the values of the independent sources in the network, we can manipulate 
the matrix equations to find the values of the current and voltage 
associated with each resistor. Remember that voltages, currents and 
resistors have analogues in other types of physical system; for example, in 
a hydraulic system we may wish to work out the pressure drops and fluid 
flow rates associated with valves that behave as hydraulic ‘resistors’. 
We can use the method described in this section for any linear physical 
system that obeys Kirchhoff-type laws, simply by interpreting the 
currents and voltages as the through and across variables of the system 
under consideration. 


The graphical approach we describe for the analysis of physical 
networks involves the construction and manipulation of large matrices — 
generally of order 2m x 2m for a network whose oriented graph has m 
edges. So, for example, when using this method to work out all the 
voltages and currents in an electrical network containing only three 
components (say a battery and two resistors), we need to handle a 6 x 6 
matrix. The solution procedure, therefore, is not particularly suited to 
hand calculation, but can be relatively easily programmed on a computer. 
We discuss a method called Gaussian elimination, which is a systematic 
procedure for solving a set of simultaneous equations in matrix form and 
which forms the basis of a number of computer solution algorithms. 


In Section 3.3, we extend our discussion to electrical networks that contain 
capacitors and inductors as well as resistors and independent sources. For 
such a network, the component equations involve not only v and i, but also 
their derivatives with respect to time. We deal with this type of 
equation by specifying certain of the voltages and currents as state 
variables. We then set up equations, called state equations, involving the 
state variables and their derivatives. We solve these equations for the 
state variables and show that the remaining network variables can be 
expressed in terms of these state variables. 


3.1 Formulation of the matrix equation 


We first consider the component matrix, and then show how to combine it 
with the fundamental cycle and cutset matrices. 


The component matrix 


Our first task is to write down the component equations of our network in 
matrix form. In Section 1 we discussed four types of component — resistors, 
capacitors, inductors and independent sources. For the moment, however, 
we restrict our attention to networks that contain only resistors and 
independent sources. Later, we discuss the modifications that must be 
made when capacitors and/or inductors are present. 
Now the component equations corresponding to the resistors are all of the 
form 

v=Ri, 
and the component equations corresponding to the independent sources 
have the form 

v=V for a voltage source, 
or 

isl for a current source. 
We begin by rewriting the component equations so that all the through 
and across variables appear on the left-hand side, and all the terms not 
involving these variables appear on the right, thus: 

lv-Ri=0 for a resistor; 

lv+0i=V fora voltage source; 

Ov+li=l for a current source. 
Note that we have included the coefficients 0 and 1 where appropriate — 
this is because our next step is to rewrite the component equations in matrix 
form. We illustrate this procedure by two examples. 


Example 3.1 
Consider the following network and its oriented graph: i 8: 
The component equations for the resistors are: 


maRi 


0, = Reis 
which we can write as 
1m, -Ryi,=0 
1v,-R,i, = 0. 
We can write these equations in matrix form as follows: 


1 0 '-R 0 ][x%]_f0 i 
0 1/0 -R ilo} [0 1) 


u) 
iy 
There are 2 resistors, so on the left we have a 2 x 4 matrix of coefficients, 
multiplied by a 4 x 1 column vector of edge variables, and on the right we 
have the 2 x 1 zero column vector, 
There is only one independent (current) source I, so we have only one 
further component equation, 
i=], 
which we can write as 
00, + 1i, =1. 


We can write this equation in matrix form as follows: 
o af] (2) 
1 


We now combine the matrix equations (3.1) and (3.2). To do this, we 
incorporate all the edge variables into a single column vector z as follows: 
We list the currents in the same order 
voltages as the corresponding voltages; at this 
stage, the order chosen is 
ze 


unimportant. 
currents 
4 


Thus we obtain the matrix equation 


We have shaded the partitions 
arising from equation (3.1). 


Example 3.2 
Consider the following network and its oriented graph: 


— oe 


The component equations for the resistors are: 
v,=R,i, or 1v,-R,i,=0 
v,=Ri, or 1v;-R;i;=0 


Ug=Rgig or 1v,—Ryi,=0. 
We can write these equations in matrix form as follows: 


1000 00 0 0}R0 00 0 0 0 0) fxm] fo 
0 10 0 0 0 0 010-R: 0 0 0 0 0 0] Ja} Jo 
0010000 0:0 0-0 00 00] {x%| Jo 
0 00100 0 0/0 0 0-R0 0 0 0 | Js} 0 
000 0121 0 0 O10 0 0 O-R 0 0 0 Mix, |=} 0} (31) 
0000010 0:0 0 0 0 0-R, 0 0] |x] }0 
00 000 01 0!0 0 0 0 0 0-R 0 | |x} 0 
0000000 1:0 0 0 0 0 0 0 -Rs| |v} [0 

py 

is 

ig 

is 

ig 

i 

ig 

ig | 


There are 8 resistors, so on the left we have an 8 x 16 matrix of coefficients, 
multiplied by a 16 x 1 column vector of edge variables, and on the right we 
have the 8 x 1 zero column vector. 


There is only one independent (voltage) source V, so we have only one 
further equation, 


v1, =-V, 
which we can write as 
10, +0i,=-V. 
We can write this equation in matrix form as follows: 
fi WF] =[-V} (3.2) 
h 


Combining the matrix equations (3.1) and (3.2) as before, we obtain the 
matrix equation 


10000000 0;-R 0 0 0 0 0 0 0 
010000000; 0 -Rk; 0 0 0 0 0 0 
001000000:0 0 -Ry 0 0 0 0 0 
qoo00100000!0 0 0 -Rs 0 0 0 0 
000010000; 0 0 0 O =-Ke 0 0 0 
900001000; 0 0 0 0 0 -Ry O 0 
000000100: 0 0 0 0 0 0 -Rg 0 
0o00000010:0 0 0 0 0 0 0 -Ry 
000000001; 0 0 0 0 0 0 0 0 


The above examples show how we deal with resistors, independent 
current sources, and independent voltage sources. Note that, in each 
example, the matrix of coefficients on the left is very simple — most of 
the elements are zero, and each non-zero element is either 1 or — R;, for 
some resistance R,. 


More generally, for a given network of this type for which the oriented 
graph has m edges, we write down all the component equations and then 
rearrange them in matrix form 


1 is the mx 2m matrix of coefficients which appear in the 
component equations; 
y.] is the 2m x 1 column vector of all the voltages, followed by 
i | all the currents, listed in the same order; 
f is an m x 1 column vector containing terms from the right- 


hand side of the component equations for independent sources 
(such as I and -V). 


In fact, most of the elements of the matrix T are zero, Such a matrix is said 
to be sparse. 


Note that the reference arrow for 


edge 1 of the oriented graph points in 
the same direction as the arrow on the 
voltage source in the circuit diagram. 
Because of the different conventions 
used in the two cases, the component 


equation has a minus sign in it. 


esoooocce 


pocecescoos 
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Problem 3.1 


Consider the following network and its oriented graph: 
ee} 
(oe 


Assuming that R,, R;, ..., R, and V are known, write down the component 
equations in matrix form. 


We now have the component equations in matrix form. These equations 
describe the constraints on the voltages and currents in the network that 
are imposed by the components only. But we know that the voltages and 
currents in’the network must also obey Kirchhoff’s laws. Our next step, 
therefore, is to combine the information contained in the component 
equations with the information about how the components are connected 
together which is contained in the fundamental cycle and cutset equations. 


The fundamental cycle and cutset matrices 
We saw in Section 2 that the fundamental cycle and cutset equations can be 
expressed in matrix form — namely, 

Cyv=0 and Dri=0, 
where C, and Dy are the fundamental cycle and cutset matrices associated 
with a given spanning tree, and v and i are column vectors of edge voltages 
and currents written in the same order as the columns of the matrix. 


For example, if we take the oriented graph of the network in Example 3.2, 
then the matrix equations are: 


1-1 1 0 131 0 0 O}fx] fo 
1 -1 0 0 110 1 0 O}}y |_lo 
0 0 -1 1 0:0 0 1 Ol/y | Io 
0 1 0 0 030 0 O 1]/ x] jo 

c i Mm 
6 
05 
V9 
Y 
T i) if 
(or v= 0 
and 
1 0 0 0 Of-1 - 0 O}fis) fo 
0 2 0 0 O11 1 0 -2I1%| Io 
0 0 1 0 01-1 0 1 Offi |Fo 
0 0 0 1 0:0 O -1 Ollig| Jo 
0 0 0 O 11-1 -1 0 O}lK! Jo 
I; D “Vig 
is 
iy 
uy 
i T + 
D; i=0 


The oriented graph for Example 3.2 is 
the same as the oriented graph for 
the example on page 24. We choose 
the same spanning tree, so these 
matrices are the matrices 
summarized on page 40. 


In general, it is sometimes convenient to write these equations in 
abbreviated form as 


[ClIn-nsa]v =0, 
where C is an (m—n +1) x (n -— 1) matrix, and 
[In |D]i =o, 


where D is an (n —1) x (m —n + 1) matrix. 


Combining the matrix equations 


We now have three matrix equations, corresponding to the component 
equations, the fundamental cycles, and the fundamental cutsets. These 
matrix equations are: 


component equation: 17] =f, 
fundamental cycle equation: C;v=0, 
fundamental cutset equation: D,i=0. 


Our aim is to solve these equations — that is, to work out the voltages v 
and currents i that satisfy simultaneously the component equations, 
Kirchhoff’s voltage law and Kirchhoff’s current law. Because we want to 
solve these equations simultaneously, we start by combining them into just 
one matrix equation of the form 


“(af 


This is possible, because we can write the left-hand side of each of the 
three equations as a matrix multiplying a vector of voltages and currents. 


The component equation is already in this form: 
T 7 =f 
The fundamental cycle equation can be written as 
[Cr 10] Be 
¢10}}--|=0, 
and, similarly, the fundamental cutset equation can be written as 


[0\D,] 


At this stage, a word of caution is necessary. The voltages which occur in 
the column vector v in the fundamental cycle equation appear in a 
particular order determined partly by the choice of spanning tree. The 
currents which occur in the column vector i in the fundamental cutset 
equation appear in the same order as the voltages in the fundamental 
cycle equation. It follows that if we want to combine these equations with 
the component equation, we must put the voltages and currents occurring in 


the column vector [¥] in the component equation into the same order as the 
voltages and currents in the other two matrix equations. 


Fortunately, this can be done very easily. All we need to do is to take the 
resistances R, and the zero which occur in the diagonal of the right-hand 
portion of T in the component equation, and write them in the required 


Note that C=B,! BJ? 
and D = B, 'B. in the notation of 
Section 2. 
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order. For example, in the resistor network of Example 3.2, the voltages 
appear in the fundamental cycle matrix equation in the order 

Dy Vay Dzy Vy Uy Voy U5y Voy Vy 
and similarly for the currents in the fundamental cutset matrix equation. 
The resistances and the zero on the diagonal in T must therefore appear in 
the corresponding order: 


Ry Ry Ry Ry Ry Ry Re Ro, 0. 
Once we have written the voltages and currents in each equation in the 
same order, we can combine the three matrix equations into one. 


We have 
ae =f, 
Cc, 10 
[evil] 
[0 ; Delf v ]=0, 


H_ is a 2m x 2m matrix containing component, cycle and cutset 
information; 


is a column vector containing all the voltages and currents; 


y isa column vector consisting mainly of zeros, but also including 
terms such as J and -V (source functions). 


The resistor network in Example 3.2 has 9 components (8 resistors and one 
independent voltage source), so the matrix H corresponding to this 
network is an 18 x 18 matrix. Writing the matrix equation for the network 
out in full, we get the following: 


100000 0 0 Oi-R 0 0 0 0 0 0 olf] fo 
0 170 00 0 0 0 OF 0-0 0 0 0 0 O|/y%| |0 
00100000 0:0 0-R 0 0 0 0 Olfv,| | 0 
0) 8) 0) 1 0) Otomo 0 0 0 -Rs 0 0 0 Ofjv| | 0 
O° GG: ‘Ol 1 10 Gs eorsot 10) <0) <0! = Re 0 0 Olly] | 0 
0 0 0 0 0 £0 0 10% o) 0) O 0 © 0 0 Ol|m%}] |o 
(0; 0 (0) 50° (0) 10) 140) Os Joo Jo) 0. 0 -Rs 0 0|| v5 0 
0! 0 0) 10) 0 Ch (0) HOOF 0 0) 0 0 0 -Ro O}|v] | 0 
Le aL eR a a EE tS) A VP He 0 .0...0)) a) _|-V 
Teal dU Oited con0 0000 0 0 Ol/i}=\0 
1-1 0013010 0.0) 00 0 0 0 O}/%] | 0 
00-2 10!001 0000 0 0 Olli] |o 
GA 208 LO ORO Ole Oe yOu 0.0.0" 0..0...0}] is) |.0. 
00000 00 0 0 0 10 O]/i%} | 0 
G0) 0! (ON 0 0 0 100 1 0 =I li} 10 
00000 eng Oued 0 1 Olfis| | 0 
0 0 6 oO 0 0 0 001 0-1 Ollis| | 0 
00000 a) 00 0 -1 0 ojl4] Lo 


We have partitioned each of the matrices Cy and Dy into two, to reveal 
the identity matrices they contain. 
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Problem 3.2 


Consider the following network and its oriented graph: 
A 4 


1G Re []s ual] }}a 


B B 


Taking the left-hand edge as the spanning tree, find the fundamental 
cycle and cutset matrices C; and D,. Use these, together with the matrix 


on page 44, to find the corresponding matrix equation Hx = y. 


Problem 3.3 


Use the fundamental cycle and cutset matrices given on page 41 to find the 
corresponding matrix equation for the network in Problem 3.1. 


3.2 Solving the matrix equation 22... 


The matrix H and the column vector y contain all the information we need 
to be able to work out the 2m voltages and currents in the vector x. One 
approach to solving the matrix equation is to rearrange the equation so 
that we express x in terms of H and y. We have 


Hx=y. 
In formulating the component, cycle and cutset matrix equations, we have 
taken care to ensure that the 2m equations in 2m unknowns that describe 
the network are linearly independent. This means that the matrix H has 
an inverse H~!, If we multiply each side of the above equation by H™, we 
get 


H"Hx = H"'y. 

Remember that the product H™H is equal to the identity matrix I,,,, 50 
1,,x=H'y 

or simply 
x=Hy, 


From a theoretical point of view, the problem is now solved — in 
principle, we can go on to calculate H™ for the network and hence work out 
each voltage and current in x. However, even if we use a computer to 
perform the inversion and subsequent calculations, this is generally not 
the most efficient approach. One reason for this is that we are interested 
in H™ only as a means to an end, What we are actually interested in are 
the voltages and currents in the network. So if we can work out the values 
of these variables without having to work out all the elements of H™, 
then we can reduce the computational effort and hence the time required 
to analyse the network. 


Instead of considering how to compute H™!, we briefly discuss a technique 
that forms the basis of many computer-based methods for solving a set of n 
linear equations in n unknowns. This technique is called Gaussian 
elimination. 


This is the network discussed earlier 


in Example 3.1. 
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Gaussian elimination 


Gaussian elimination is a systematic way of solving a linearly 
independent system of n simultaneous linear equations in n variables by 
eliminating the variables one at a time. This is done by performing certain 
allowable operations on the rows of the matrix of coefficients so as to 
reduce all the elements below the main diagonal (the diagonal from top 
left to bottom right) to zero. We can then successively solve the equations 
one at a time, starting with the last equation, until we have all the 


required variables. 


The operations we allow on the rows of the matrix are: 


(a) interchange two rows; 


(b) multiply any row by a non-zero number; 


(c) add any row to any other row. 


If we also carry out these row operations on the column vector y, then each 
of them corresponds to an operation on the corresponding equations which 
leaves the solution of the equations unchanged. For example, operation (a) 
corresponds to interchanging two equations, which clearly does not affect 


the solution. 


Example 3.3 

Consider the set of equations 
2x, +3x,- x4=5 
4x, + 4x -3x, =3 
2x, -3x, + x, =1. 


These can be written in matrix form: 


2 3 -Iifm] [5 
4 4 -3//x,/=/3 
2-3 Ijixs 1 


We want to work out the values of x,, x, and x, without inverting the 3.x 3 
matrix. We do this by performing the allowable row operations with the 
aim of changing the values of all the elements of the matrix below the 


main diagonal to zero. 


The result is a matrix equation of the form 


yy Sia Sis ah 
0 ayy 35!) 2 |=| Yo 
35 |1%3| LYs 


where a), # 0, ay) #0, ay, #0. 


This is equivalent to the three equations 


MX + 2% + Ags =Y, 
Daghr + Unxks = Yo 
3X3 = Ya 
From the last equation, 
5 =Y3 (a3 #0) 


we obtain 


From the second equation, 

MypXy + AygXs= Yo (tgp #0) 
we obtain 

% tes 


I 


Note that operations (b) and (c) mean 
that we can add (or subtract) any 
multiple of any row to (or from) any 
other. 


Our aim is to describe the method 
of Gaussian elimination as simply 
and as briefly as possible. The 
example used is inevitably rather 
artificial. In practice, if solving this 
system of equations by hand, we 
would, of course, proceed differently, 
Moreover, when solving a system of 
equations by computer, we would 
design the program so that rounding 
errors are minimized. 


It is possible to do this, and to obtain 
non-zero elements on the diagonal, 
because the equations are linearly 
independent. 


Substituting x, =, /aj,, we get 
oa "Ns 
842 Ap2835 
Similarly, we can work out x, by substituting for x, and x, in the first 
equation, 


This method of solution is called back-substitution and can be extended 
to systems of linearly independent linear equations of any order. 


Let us consider how we use Gaussian elimination to reduce the 
3x3 matrix of our example to triangular form. The matrix is: 


2 3 -1 
4 4-3 
23 1 


We start with the first column. We require that all the elements in this 
column except the top left be zero. 
If we multiply the first row by — 2, we get 

=4 6 2 


Adding this to the second row reduces that row to 
0 -2-1. 

Turning our attention to the third row, we subtract the first row to give 
0-6 2 


5 
We must also remember to carry out the same operations on the vector | 
1 


on the right-hand side of the matrix equation. Thus the first element 5 of 
this vector is unaffected, the second element 3 becomes 


3+(-2x5)=-7, 

and the third element 1 becomes 
1-5=-4, 

The matrix equation can now be written as 


2 3 -1)fmn 5 
0 -2 -1|)x2.|=|-7 
0 -6 2}/x3}] |-4 


Our operations to modify the first column using the first row are complete. 
We now consider the second column and the third row. Multiplying the 
second row by —3 and adding the result to the third row, we reduce the 
third row to 


O10) 


so the second column has the required form. 


The corresponding term ~ 4 on the right-hand side becomes 
-4+(-3x-7)=17. 
Thus the final matrix equation is 


2 3 -1\/x 5 
0 -2 -1]} x, |=|-7 
0 0 5]ix5 17 


From the third equation, we find that 5x, = 17, or x, = 17/5. Substituting 
this value for x, into the middle equation and solving for x,, we get x, =9/5. 
Finally, knowing x, and x,, we can solve the first equation for x,, and obtain 
x, =3/2. 


This is to eliminate the term in x, 
from the second and third equations. 
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Using Gaussian elimination to convert the matrix to triangular form, and 
back-substitution, we have solved systematically the system of linear 


equations without having to invert the 3 x 3 matrix. 


We now apply this technique to the solution to a network problem. 


Worked problem 


Find the currents and voltages in each part of the following network: 


Solution 
The matrix equation, obtained in the solution to Problem 3.2, is 
fo Oo Of1 © 0 )fu) [F 
ee) ' O -R, 0 |/m] |0 
0 0 110 0 =Rs|}v3|_]0} 
-1 1 0:0 0 O jl] "Jo 
1 0 1;0 0 Ol 0 
0; Oya a ilsil (0 
Interchanging row 1 and row 4, and row 5 and row 6, we get 
ft oT Hor 30; On 0 
0D 26 ' 0 -R, 0 || 0 
0 0 110 0 -Rs|}v3|_]o} 
0.0 O8L 0 ONT 
5A me TNE fay Os es 0 OF 0 
-1 0 110 0 0 jfis} [o} 
Subtracting rows 1 and 3 from row 6, and row 4 from row 5, we get 
[=l 2 OHfG 10) 0 ] »] [0 
Cee te) \ 0 -R, 0 {|x 0 
0 0 110 O -Rs/}v3|_}0 
OO ONT DO Od 
O 20) ONO, AL ais -I 
Lo -1 0:0 0 Rs Jlis| Lo 


Finally, after adding row 2 and R, x (row 5) to row 6, we get the matrix 


equation in the required form: 


Successively solving these equations, beginning with the last one, we 


obtain 


a 


ecoose 


1 
1 
0 
0 
0 
0 


(Ry +R3)is =-RyL, 


i, +i,=-I, 
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ocoroo 


8888 8 8 


corooco 


0 0 Dy 
-R, 0 0 

LS a ie 
0 Oo Ha)” 
i 1 iy 

O R2+Rs]| is 


i, =—RyI/(R, +R), 
i, =—Rgl/(R, +R), 
i,=1, 

D3 =—RgRgl/(R, + Ry), 
0, =—RyRgl/ (Ry + Rs), 
2, =-R,RyI/(R, +R,)- 


0 
0 
0 
i 


ih 
Ral 


Since i, and i, are negative, the 
current must flow from B to A in this 
part of the network. 


We have now achieved our goal. Given a network containing resistors and 
independent sources, we can work out the voltage and current associated 
with each edge of the oriented graph of the network. 


Remark We have formulated the matrix equations so that they can be 
solved straightforwardly for cases in which we know the values of the 
source voltages or currents and wish to calculate the voltages and currents 
in all the other components. But suppose that we want to solve the 
converse problem and work out, say, the value of a source voltage which 
will result in a specified current flowing in one of the resistors of the 
network. To do this, we would have to rearrange the matrix equations, 
Essentially, we would have to treat the source as if it were a resistor and 
treat the resistor which is to have a specified current flowing in it as if it 
were a source. We mention this because such problems can arise in 
electrical network design, but we shall not discuss such a case in detail 
here. 


All the examples we have discussed in this section have involved only 2- 
terminal components. However, the methods of solution we have given can 
be applied in a similar way to systems involving components with more 
than two terminals — for example, transformers or transistors, The cutset 
and cycle equation matrices are of the same form, but the component 
equation matrix is different. 


mg be inhaled ong 
: Hy Ascsat (Mob equGor Har at KA my onl 
3.3 State equations trey be ciel to spain Te Sls. ich.) 

Up to now, we have confined our discussion to networks that contain only 

resistors and independent sources. For such networks, the component 

equations are particularly simple, being proportionality relations of the 

form v = Ri. 


We now turn our attention to networks involving capacitors and inductors. 
For such networks, the component equations involve not only v and i, but 
also their derivatives. Our aim is to show how the matrix equations can 
be adapted to take the derivatives into account. 


In Section 1 we saw that the component equations of a capacitor and an 
inductor are 


._ ado di 
i=CG, and v=La, 
respectively, where C is the capacitance and L is the inductance of the 


component, Using a dot to denote a derivative with respect to time t, we 
may rewrite these equations in the form 


i=di and isto. 
We shall use these equations to obtain the matrix equation in the form 


Hx=y 


where 


x is a column vector involving all the voltages and currents; 


x is a column vector involving their derivatives; 
is the ‘source function vector’ we introduced earlier; 


x 
H_ isa large matrix involving the components and the fundamental 
cycles and cutsets; 


K_ is a matrix isolating the derivatives of the edge variables. 


We illustrate the ideas involved by considering the following simple 
network and its oriented graph: 
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A 
1 8: ile 3 


i L q 3 
B B 
The component equations for this network are 

i = I(t), 

nye 

=F ily 

eee! 

1h =[ 05. 


We write these equations in matrix form as follows: 


0 0 1 0 {0 ojfm] fm] fo 0 0 Oo; 0 ofa 
0 WL 0 0 {0 Ollys|=|0|+/0 0 0 0:0 1I\% 
0 0 0 10:0 O}]i) LO 0 0 0 0:1 Oh 
By B 
"a 2 
iz) i) 


Notice that the variables whose derivatives occur in the component 
equations (v, and i,) appear last of all in the column vectors x and x. These 
variables are called state variables. The values of the state variables 
represent the state of the system. 


We now wish to combine this equation with the fundamental cycle and 
cutset equations for the network, These are the same as the fundamental 
cycle and cutset equations for the network of Example 3.1, and can be 
written as follows: 


The shaded columns correspond to the state variables v, and i,. 


Rearranging the columns so as to put the edge variables into the same 
order as in the component equation, we get 


1 0 0 0 v] fo 
Sea) 0 »| =|0 
oe 1)| i,| [0 


OO; 
Combining this with the component matrix equation above, we obtain the 
required equation: 


0 0 1 0 {0 O}fy] fi] fo 0 0 O10 O}fy 
0 YL 0 0 10 Ol}o} | 0] Jo 0 0 O10 1)],, 
0 0 0 YC!0 Oj} |_} 0 |,/0 0 0 O}1 OF], 
-1 0 0 011 Offi |"| 0 |*}o 0 0 Oo off” 
-1 1 0 0 10 Of/%| | 0] jo 0 0 Of0 Offa 
0 © 1 2 10 INK 0} [0 00 010 Ofly 

t 7 t ‘t 7 

H x = y + K x 


‘This matrix is the lower half of the 
matrix in the solution to 
Problem 3.2(a). 


Our next step is to manipulate these matrices so that we can write down 
two state equations for the network. The state equations are first-order 
differential equations that relate the derivatives of the state variables 
to the state variables and the independent sources. It can be shown that 
the two state equations are sufficient to describe completely the dynamic 
behaviour of this network. In other words, any voltage or current in the 
network may be expressed in terms of the state variables v, andi,, their 
derivatives 3, and é, and the independent current source I(t). 

To solve the matrix equations for v, and i;, we use Gaussian elimination. 
But, because we are interested only in v, and i,, we need concentrate only on 
the rows of the H matrix which involve these variables, that is, on rows 4 
and 6. 


We want the final matrix equation to be of the form 


ry cy 
U3 az ? ? Uo) 
ale h 
ig)7|...)*]0...0..0.. ip 


Those regions denoted by ? are to be determined. 


We wish first to eliminate the term —1 in the first column of row 4 of H. 
We do this in two stages. First we subtract row 5 from row 4, giving the 
new row 4 of H: 

O-1 001 0. 
A new term, -1, has now appeared in the second column, but we can 
eliminate that by multiplying row 2 by L, and then adding it to the new 
row 4, giving 

oO 0) Ol 10 F 10: 
Next, we perform the same row operations on the right-hand side, and the 
matrix equation becomes: 


i¢t)| [0 0 0 0} [o, 
0| |o o o 1|| ds 
0 |,/0 0 0 OnE 
eR Vl eee Alle 
07) jo oo 0} |e, 
Of 102 Sa TO! (OF OF (Oye 


The first four numbers in row 4 are now zero, as required, and we turn our 
attention to row 6. 


Adding -1 x (row 1) and —C x (row 3) to row 6, we get 


m] [1] fo 0 oy 
03 0 0 0 Ly 
|_| 0 0 0 n 
4] |0.]"]o.0 0. h 
WY 0 0 0 d 
| |-1()} [0 0 I 


This completes the elimination procedure. From rows 4 and 6 of the matrix 
equation we can write down the two state equations relating the two state 
variables of the network: 


= Lis 


i5=-I(t)- Coy 


Notice that row 4 has a 1 in column 5, 
corresponding to the 5th element in x, 


namely v,, and that row 6 has a 1 in 
column 6, corresponding to i,. 
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or, in the more usual form, 
Feo! 
B= % 


b=-Zi, = z I(t). 
The general matrix form of the state equations is 
x=Ax+u, 
where 
A isakxk matrix; 
x is a column vector of the k state variables of the system; 
x is a column vector of their derivatives; 
u__ isa column vector involving the independent sources. 


For example, the state equations in our example can be written as 


|_| 0 -/C]/ 09] , |-(UC)I(E) 

4] [WL 0 | lis 0 
We now solve the state equations. The exact form of the solution depends 
on how the current I(t) varies with time, on the initial conditions of the 


problem (that is, the values of v, and i, at time t = 0), and on the values of 
Land. 


In this case, the solution has the form 
v,(t) = A cos wt + B sin wt + fit), 
i,(t) = A, cos wt + B, sin wt + g(t), 


where A, A,, B and B, are constants, @ = (LC)/* and f and g are functions 
which depend on I(t) and its derivative. 


In fact, 
A, =-BCo, 
B, = AC, 


git) =- Cf’ (t) - (8), 
and f is a function satisfying the differential equation 
1 1s 
f'O+ FEAN+ EMO =0. 


Finding the other edge variables i,,v,, i, and v, is now straightforward, 
since we can obtain them in terms of the state variables v, and i, and their 
derivatives. To do this, we look back at the matrix equation (*). The first 
four equations are 


i, = I(t), 
ts 
L378 
etre 
Cary 
-v, + Ly =0, 


and i,, 03, i, and v, can easily be calculated from these equations. 


Problem 3.4 


(a) Express the following component equations in matrix form: 
y do. di. i 
=k i =Cap, Vy=Lapy 4 = Rig. 


(b) Write down the H-matrix equation for a system which has the 
component equations given in part (a), and the following 
fundamental cycle and cutset equations; 


1 11 Olf~]_fo 10 -1 0}[4]_fo 
F 10 1 2% -(5| by b 1-1 | 3 -(1] 
2 i 
04 ig 


(c) Find the state equations, 


The procedure carried out above for our simple example can be extended to 
any network made up of capacitors, inductors, resistors and independent 
sources, or their analogues in other types of physical system, The number 
of state variables, and hence state equations, needed to describe the 
dynamic behaviour of a system completely is related to the number of 
energy storage components in the system. In our simple example, we had 
two such components: a capacitor, which stores energy in the form of an 
electric field, and an inductor, which stores energy in the form of a 
magnetic field. We were able to write down the general form of the 
solution to the state equations, but for more complex networks the state 
equations would be solved numerically using a computer, There are various 
computational techniques available to tackle the problem, but a 
discussion of their relative advantages and drawbacks is beyond the scope 
of this unit. 


3.4 Computer activities 


The computer activities for this section are described in the Computer 
Activities Booklet. 
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Further reading 


Further information about electrical networks may be found in: 
C. A. Desoer and Ernest S. Kuh, Basic Circuit Theory, McGraw-Hill, 1969; 


S. W, Director, Circuit Theory: A Computational Approach, John Wiley and 
Sons, 1975; 


W. A. Blackwell, Mathematical Modelling of Physical Networks, 
MacMillan, 1968. 


Further information about physical networks in general may be found in: 
P. H. Roe, Networks and Systems, Addison-Wesley, 1966. 

A useful account of matrix theory, illustrated with many examples, is 
given in: 


F, Ayres, Schaum's Outline of Theory and Problems of Matrices, McGraw- 
Hill, 1962. 
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Exercises 


Section 1 


Components, and through and across variables 

1,1 Classify each of the following as a through variable or an across 
variable: 

(a) the fluid flow rate for a hydraulic valve; 

(b) the force acting on a mass; 

(c) the derivative of the voltage for a capacitor. 

1.2 Identify the through and across variables, and write down the 
component equation, for each of the following components: 

(a) an electrical inductor; 

(b) aspring; 

(c) a transformer (considered as a 3-terminal component). 


1.3 Consider the following system: 


(a) Draw an oriented graph representation of the system. 
(b) Apply Kirchhoff’s voltage law to each cycle of the oriented graph. 
(c) Apply Kirchhoff’s current law to each vertex of the oriented graph. 


1,4 Repeat Exercise 1.3 for the following system: 
4 B 


D 


1.5 Repeat Exercise 1.3 for the following system: 


Graphical representation 


1.6 Decide whether each of the following statements is TRUE or FALSE. 

(a) Inan oriented graph, a 6-terminal component can always be 
represented by 6 vertices and 5 edges. 

(b) Inan oriented graph representation of an electrical network, the 
arrow next to an edge always points in the same direction as the 
actual direction of the current flow in the corresponding component. 

(c) Inanoriented graph, a 5-terminal component can be represented in 
exactly 5 different ways (ignoring edge directions). 


1.7 The following graphs can both be used to represent the same 


4-terminal component: 
a1 8B 4 B 
Ve alt He 
> ae a Dia se 
(a) Express each of the voltage variables v,, v, and v, in terms of v,, 0, 
and 23. 


(b) Express each of the current variables i,, i, and i, in terms of i,, i, and i,. 


1.8 Draw an oriented graph representation of a 5-terminal component. 
How many component equations are associated with such a component? 


1.9 Draw an oriented graph which represents the following electrical 
network: 


(a) Construct the dual network. 

(b) Write down the component equations of the dual network, given that 
the soronet ations a origiagl network are 

; di 

v,=3, % =i, vg Sig 1,= 47, gn 6g? 

F “x Nee 

Section 2 

Fundamental cycles and cutsets 


2.1 Consider the following oriented graph, in which a spanning tree is 
indicated by thick edges: 


(a) List the fundamental cycles and write down the corresponding 
voltage law equations. 

(b) List the fundamental cutsets and write down the corresponding 
current law equations. 

(c) Find the fundamental cycle matrix C,. 

(d) Find the fundamental cutset matrix D;,. 


2.2 Anoriented graph has 8 vertices and 12 edges. Write down: 
(a) the number of edges in a spanning tree; 

(b) the number of edges in each co-tree; 

(c) the number of fundamental cycles; 

(d) the number of fundamental cutsets. 
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Obtaining the fundamental cycle and cutset equations 


2.3 Consider the oriented graph of Exercise 2.1, with vertex A as 
reference vertex. 


(a) Write down the reduced incidence matrix, and partition it into a tree 
part B, and a co-tree part B.. 


(b) Use the method described in Section 2.5 to find the inverse of the 
matrix By. 


(c) Use the matrices By and B. to find the fundamental cycle matrix C, 
and the fundamental cutset matrix Dy. i 


(d) Calculate the matrix product C,D/. 
2.4 For each of the following oriented graphs, write down the 


fundamental cycle equations and the fundamental cutset equations, using 
the spanning tree indicated by thick edges. 


Hint These oriented graphs are the duals of those on pages 40 and 41. 
Use the results stated in Solution 2.10. 


Section 3 


Formulation of the matrix equations 


3.1 The following diagram shows an electrical network and an 
associated oriented graph: 


a B The resistances are given in ohms (Q). 


6vt 4a 


©; 


(a) Using the spanning tree shown, find the component equations and the 
fundamental cycle and cutset equations. 


(b) Solve the equations to find all the currents and voltages in the 
network. 


3.2 
(a) Express the following component equations in matrix form: 
B= Ri, V2 =3, Vy=Rgig, % =Ryig- 


(b) Write down the H-matrix equation for a system which has the 
component equations given in part (a), and the following 
fundamental cycle and cutset equations: 

1 0 -1 O}fi]_fo 
0 1-1 Ils] [0 


1 1 1 O/fm)_fo 
0 -1 0 Ilo3|=|0 
02 Qn 
04. ig 
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Solving the matrix equations 
3.3 Use Gaussian elimination to solve the following matrix equation: 
100 1fx 3 


120 Offx|_| 3 
16 1 illx|| 0 
00 3 -Af{x4] [-13 


State equations 


3.4 The following diagram shows an electrical network and an 
associated oriented graph: 


Using the spanning tree shown, find: 
(a) the fundamental cycle matrix; 
(b) the fundamental cutset matrix; 
(c) the component equations; 

(d) the H-matrix equation; 

(e) one of the state equations. 
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Solutions to the exercises 


Ast 
(a) 


1.2 
(a) 


(b) 


(c) 


13 


(a) 


(b) 


(c) 


1.4 


(a) 


(b) 


through variable; (b) through variable; (c) across variable. 


The through variable is the current i and the across variable is the 
voltage v. The component equation is 
di 
v=Lap 
where L is the inductance. 


The through variable is the applied force f, and the across variable 
is the extension x. The component equation is 


fakx, 
where k is a constant. 


Since the transformer has three terminals, it has two components. 
The through variables are the currents i, and i, flowing in the input 
and output windings, and the across variables are the voltages v, and 
v, across these windings. The component equations are 


v, =(n,/n,v, and i, =-(n,/M,)in, 
where n, and n, are the numbers of turns in the two windings, 


2 
4 B 3 
. D 


system oriented graph 


cycle ABDA v, -% + v, =0; 
cycle BCDB v, ~v,+ v, =0; 
cycle ABCDA v,+v,-v, +,=0. 


vertex A -i, +i, =0; 
vertex B  -i, —i, +1,=0; 
vertex C -i, —i, = 0; 
vertex D i, +i, +i; =0. 


2 (4 
t Is y 34 
5 
D eae) 
system oriented graph 


cycle ABDA v,-v, -v, = 0; 
cycle CDC “Us +0, =0. 


vertex A i, +i, =0; 
vertex B -i, -i, =0; 
vertex C =i, —i; = 0; 
vertex D -i, +1, +i, +i; =0. 


system. oriented graph 


(b) cycle ABA (left) 0, +0, =0; 
cycle ABA (right) -v, +, =0; 
cycle ABA (outer) -v, + 0,=0. 


(c) vertex A -i,-i, -i, =0; 
vertex Bi, +i, +1, =0. 
1.6 
(a) TRUE; (b) FALSE; (c) FALSE. 


1.7 
(a) Clearly v, =v; and v,=—v,. 
The potential difference across the terminals D and A is 
Renee ect tiae a 


(b) Applying Kirchhoff’s current law to the vertices A, B and C, we 


have 
i,=-i, i,=i,-i, ip+i,=i,-i, 
giving 
ige-iy ipeig-i, igi, -i. 
1.8 x 
There are several possibilities — for example: ~2\ Y, 
All oriented graph representations are trees with four edges, so there are T 


four component equations associated with a 5-terminal component. 


1.9 


oriented graph dual graph dual network 


Notice how the cycles/cutsets in the dual oriented graph correspond 
to the cutsets/cycles in the original oriented graph; for example, 
look at the edges marked 1, 5, 4 and those marked 1, 5, 2 in each 


oriented graph. 
= © GS © 


; di 
sie tesa lice sett 
a. Sos 


iou, 


fundamental fundamental 
chord _ cycle cycle equation 
BE BEDCB B, +0, -1-4,=0 
CB CBC fl 2, +0,=0 

¢c 

CE CEDC U5 +0, —0, =0 
AE AEDCBA 04 +0, - 1 - y+ 0,=0 
(b) 

fundamental fundamental 
branch _cutset cutset equation 

ay, 

BA (BA, AE} i, ig =0 
BC (BC, CB, BE, AE} igi +i tig =0 
cD (CD, CE, BE, AE} igtis +i +i=0 
ED (ED, CE, BE, AE} i,-i5—i,-ig =0 
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BC 
-1 


2.2 
(a) 7%; (b) 5; 


23 
(a) The reduced incidence matrix is 
BA BC CD ED BE CB CE 
110 O}1 10 
West (iit gl th 
a @ 2 -lo 0 0 


000 1!t0-1 
By Bo 


& 
0 
0 
) 


moos 


=1 
(b) 


Thus the matrix Bj is 


BiniGi B) 
pat l 1 1 1 
sc]O -1 -1 -1 
co}0 0 +l +1 
eojo0 0 0 1 
(c) The matrix product B;'B. is 
111 1f1-1 0 
0-1-1 -1}] 0 1 1 
000 
-1 0-1-1 


The fundamental cycle matrix is 
0 -1 
i 
0 -1 
Pt 

The fundamental cutset matrix is 
1 


Bs 0 
Ce=[-(B;'B.)' MIg]=] 5 


Dy =[14 BiB 


coro 
oroo 


0 
0 
0 
(d) CD =0. 
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Tracing a path in the tree from each vertex to the reference vertex, 
we find that the edges ED and BA are traversed in the direction of 
their arrows, and CD and BC are traversed in the opposite direction. 


coro 


(c) 5 (one for each chord); (d) 7 (one for each branch). 


Note that these answers agree with 
those given in the solution to 
Exercise 2.1. 


oroo 
rooo 


2.4 


Using the results on pages 40 and 41 and the results about dual oriented 
graphs, we can immediately write down the equations as follows. 


For the fundamental cycle equations, we replace i by v in the fundamental 
cutset equations and obtain the following: 


graph (a) graph (b) 

U3; -0; —U, =0 U,-0, =0 

0, +, +0; -0, =0 0; +% -v, =0 
U,—U, +05 =0 U, +0 —v, =0 

0, —% =0 U5 +, -v, =0 


U4 -¥, —U; =0 


For the fundamental cutset equations, we replace v by iin the fundamental 
cycle equations and obtain the following: 


graph (a) graph (b) 
i tigti, +i,-i, =0 i,-i,-1,=0 
i, +i, +i, -i, =0 is +i,-i,=0 
+i,-i,=0 i, +i, +i, +i, =0 
i, +i, =0 
34 
A B 
6v{ 40 


Cc 


(a) The component equations are: 


component 1 (voltage source) v, =-6 (lL) Note that the component equation for 

component 2 (2 Q resistor) 0% =2i, (2) the voltage source has a minus sign in 

component 3 (4 Q resistor) 03 =4i, (3) it. This is because the reference arrow 

component 4 (4 Q resistor) 4 =4i, (C)) for the edge corresponding to the 
voltage source defines a positive 

The fundamental cycle equations are: voltage as one which is in the opposite 
direction to that of the voltage source 

chord cycle cycle equation in the network diagram. 

AB ABCA v, +0; +0, =0 (5) 

BC BCB 1, +0, =0 (6) 


The fundamental cutset equations are: 


chord —_cutset cutset equation 

CA (CA, AB} i, -i,=0 (7) 

BC (BC, AB, i,-i,-i,=0 (8) 
CB} 


(b) We have eight equations involving eight variables. We solve these 
equations as follows. Substituting equations (1), (2), (3) and (4) into 
equations (5) and (6), we get 

2i, +4i,-6=0, so i, =3-2i, (= i, from equation (7)) 
and 

4i,+4i,=0, so i, =-i5. 
Substituting into equation (8), we get 


i,-3+2i,+i,=0, so is=pamps. 
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Hence 


eee ns ae lag: 
i, ==7 amps; i,=gamps; i, =-Zamps; 
from equations (1), (2), (3) and (4), we now get 
v, =-6 volts; v,=3volts; v,=3-volts; v, =—3 volts. 


3.2 
(a) The component equation matrix is 
100 0/-Rk, 0 0 Ojfm] fo 
0100! 0 -R; 0 Offv,|=l0 
0010! 0 O -R,y O}/m%} jo 
0001!0 0 0 Ofjm] [3 
iy 
is 
iy 
i 


(b) Before forming the H-matrix we first rewrite the matrix of part (a) 
so that the voltages and currents in the column vector appear in the 
same order as that of the fundamental cycle and cutset matrices. The 
component equation matrix thus becomes: 


10 0 OFF © 0 © 
0 1 0 0/0 -R 0 0 
0 0 f O10 @ 0 0 
[0h O) TOF D470) 10) 70) r=Ry 
The H-matrix equation Hx = y is therefore 
fi 0 O OfS% @ OO O Tia] fo 
H 
0 1 #0 Of 0 -Rs 0 0 |fn5} Jo 
OO. 040 0) DO Ii a 
0 0 0 140 0 0 =R|}r4|=}0 
' i | JO 
\ ig] |0 
o 0 °O0 OF 1 O + O hig [0 
Cae el ol SS ST ATE li 
3.3 
We are given the matrix equation 
100 tx 3 
120 Olle! | 9 
101 tlxz/"| 0 
0 0 3 -I{I/x,| |-13 


We first modify the first column by using the first row. To do this, we 
subtract the first row from the second and third rows to give 


row2: 0 2 0) -1 
row3: 00 1 0 


The matrix equation can now be written as: 
100 1)[x, 3 


020 -1/Ilx,|_| 0 
001 offx;|=| -3 
00 3 -1|/x,] |-13 


The second column is already in the required form, so we turn our attention 
to the third column. To make the bottom element (3) zero, we must 
multiply the third row by 3 and subtract it from the fourth row. The 
matrix equation then becomes: 
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100 1)[fx 3 
0 2 0 -1/\x, 0 
0 0 1 Of}x3|~}-3 
0 0 0 -1}\x,} |-+4 


We now find the values of the variables by back-substitution. 

From the fourth equation, x, = 4. 

From the third equation, x, =-3. 

The second equation is 2x, -x, = 0, giving x, = 2. 

The first equation is x, +x, = 3, giving x, =—1. 

So the solution of the matrix equation is x, =-1, x, =2,x,=-3and x, =4. 


3.4 


chord cycle cycle equation 
CA CABDC 0, +0; +05 +0, =0 
BC BCDB U3; —U, —v, =0 
AD ADBA 14 -05-0,=0 


The fundamental cycle matrix is 
2 ¢ tie a % 
AB BD DC CA BC AD 


Cad: tet 1510 0) 
Bc} 0 -1 -1 0 1 0 


AD[-1 -1 0 0 0 1 
(b) 
branch cutset cutset equation 
AB (AB, AD, CA} mi, +i, +i, =0 
BD (BD, AD, CA, BC} -i, +i, +i, +i, =0 
Dc (DC, CA, BC} ~i +i, +i, =0 


The fundamental cutset matrix is 
AB BD DC CA BC AD 
4B} 1 0 0-1 O 1 
BoD} O 10-1 1 1 
0) 1Or or =f 10) 


(c) The component equations are: 


voltage source --V 

resistors %=Ryb, 05=Ryiz, vs =Reis 
do, 

capacitor iy =Cat 
di, 

inductor 1, =Lap 
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(d) The H-matrix equation is Hx = y + Kx which can be written as 
follows. (Note that the variables v, and i,, whose derivatives occur 
in the component equations, appear last.) 


< 
D 
S 


Scocorcoocon 
ecotorocecto 
csoooroocorncse 
eorocococoooko 
oroscoocoRSOSS 
eoeorescocoo0a 
Hrososcoscoo90S 

poe 

" 

ecoosccoco0d 

+ 
eecocecccoocca 
ceocooscsceoscoo0ne 
ecoooosoHrooccle 
ecooocoxrocoscos 


(e) We can eliminate the (-1)s in row 9 by taking 
row 9 + row 7 — row 1—(L x row 6), giving 


000000000010 
Applying these row operations to the right-hand side of the 
H-matrix equation, we obtain the state equation 


= V-Lig. 


Solutions to the problems — 


Solution 1.1 


There are several possible oriented graphs, depending on how we label 
the components and which direction we assign to each edge. One 
possibility is shown in the margin. 


Solution 1.2 
Differentiating the given equation, we obtain 


In this case, f is the through variable and u is the across variable. This 
equation is of the same form as the component equation for an inductor, 


di 
v=Li 


so the spring is analogous to an inductor, and the constant 1/k (the 
reciprocal of the spring stiffness) corresponds to the inductance L. The 
analogous quantities for the two components are shown in the following 
table. 


mechanical electrical 

component spring inductor 
through variable f i 
across variable u v 

component equation aldf opti 

ie "kat onle 


Solution 1.3 


(a) ry A 
TaN d 
B c B = (4 


Applying Kirchhoff’s current law to the vertex A, we obtain 
current flowing into A = i,. 
Applying Kirchhoff’s current law to the vertex A for the original 
graphical representation, we obtain 
current flowing into A = i, +i,. 
Hence 
i,=i,+i,. 
Applying Kirchhoff’s current law to the vertex C, we obtain 
current flowing into C = -i,,. 
Applying Kirchhoff’s current law to the vertex C for the graphical 
representation in the text, we obtain 
current flowing into C = -i,. 
Hence 
ip=i. 
The potential difference (or voltage) across terminals A and B is just 
v,: (The direction of the arrow indicates that this is considered 
positive if A is at a higher potential than B.) Looking back at the 
graphical representation in the text, we see that the potential 
difference across terminals A and B is v,. Hence 
0, =. 
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The potential difference across terminals B and C is v,. Looking at 


the graphical representation in the text, we see that the potential 
difference across terminals B and C is v, —v, (that is, the potential at 


B minus the potential at C). Hence 


Uy =U —0,. 
(b) 
7 X \ 
B i 7 
Using the same method as in part (a), we obtain: 
i,=i, +i, ip=-Ay, Vg=Vy, y=, -Yy. 
Solution 1.4 


(a) To find a graphical representation of this system, we replace each 
2-terminal component by an edge, each 3-terminal component by a 
tree with two edges, and the 4-terminal component by a tree with 
three edges. Since each tree can be chosen in several different ways, 
there are many possible ways of representing this system. One 
possible way is the following. The tree edges are shown by thick 
lines. 


component edges 
ACE (3-terminal) AC, CE 
BDG (3-terminal) BD, BG 


A 
SS : 
a CDFG (4-terminal) CE, DF, FG 
E el 
F G 


(b) For any component, the number of component equations is equal to the 
number of edges in its graphical representation. So, to find the 
number of component equations associated with the system, we 
simply count the number of edges in the graph. There are 12 edges, so 
there are 12 component equations. 


Solution 1.5 
(a) We construct the dual oriented graph as follows: 


oriented graph dual oriented graph 


We obtain the dual network shown below by replacing 


component 1 by a current source; 
components 2 and 7 by inductors; 
components 3, 4, 5 and 9 by capacitors. 
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2 
eile 
2A || 


Y vA 
network dual network 


(b) We construct the dual oriented graph as follows: 


Ss. 
ie 


oriented graph dual oriented graph 


We obtain the dual network shown below by replacing 
component 1 by a current source; 
components 3 and 6 by capacitors; 
component 4 by an inductor. 


network dual network 
Solution 2.1 
(a) cycle voltage equation 
ABCEA v, +0, +0;,+0,=0 (1) 
BDCB %-0,-0,=0 (2) 
CDEC U+0,-v,=0 (3) 
ABDCEA +4 -%, +0,+0,=0 (4) 
ABCDEA Wty +% +, +0, =0 (5) 
BDECB 4% +0,-0;-2,=0 ©) 
ABDEA vy +y +0, +0, =0 (7) 


(b) The maximum number of linearly independent equations is 3. 

For example, equations (1), (2) and (3) are linearly independent, 
since each contains at least one term not included in either of the 
other two, and so none can depend on the other two. However, the 
remaining equations (4)-(7) all depend on the first three, as follows: 

equation (4) = equation (1) + equation (2) 

equation (5) = equation (1) + equation (3) 

equation (6) = equation (2) + equation (3) 

equation (7) = equation (1) + equation (2) + equation (3). 
It follows that, if we wish to study the network, we need choose only 
equations (1), (2) and (3) (or any other 3 linearly independent 
equations), since the remaining equations give us no futher 
information. 
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Solution 2.2 
(a) There are many possible spanning trees — for example: 


A B A ¢ B Cc 
bi E D F E D F E D 
branches chords branches chords branches chords 
AB AF AB AF BC AB 
BC BE BE BC BE AF 
cD CE CE cD cD CE 
DE FA DE FA EF DE 
EF EF FA 


(b) We know that every tree with n vertices has n — 1 edges. Thus there 
are n—1 branches, and so there are m —(n-1)=m-—n +1 chords 
associated with each spanning tree. 


(c) If we add a chord to a spanning tree, we obtain a single cycle. For 
example, if we add each chord separately to the spanning tree in our 
example, we get the following cycles: 


F E F E E D F 


chord BC chord BE chord CD chord FA 
Solution 2.3 
(a) 
fundamental fundamental 
chord cycle diagram cycle equation 
B 
3 
BD BDCB 4 %4-%-v, =0 
6] 
D 


DE CDEC U, +0, —-v; =0 


EA ABCEA +0; +05 +0, =0 


The cycle ABCDEA contains the chords DE and EA. We.obtain the 
voltage law equation for this cycle by combining the cycle equations 
for these two chords, giving 


v, +0, +0,+0, +0, =0. 
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(b) 


fundamental fundamental 


chord cycle diagram cycle equation 
Ac 3B 
y 
AB ABDEA ‘ 4 % +y +0, +0, =0 
7 
D 
3 
cD BCDB 4 2, +0,-0, =0 
6 
D 


CE BCEDB 0; +05 -0,-0, =0 


The cycle ABCDEA contains the chords AB and CD. We obtain the 
voltage law equation for this cycle by combining the cycle equations 
for these two chords, giving 


U, +0, +V, +0, +v, =0. 


Solution 2.4 


(a) The chords are BD, DE and EA, and the fundamental cycles are 
BDCB, CDEC and ABCEA. The fundamental cycle matrix is 
therefore 


AB BC CD CE BD DE EA 


Bo} 0) =) <1 (0 1 0) “oO 
DEMO: 0) T= 1 100 ay 8 
EA 0 A oF eg 


and the fundamental cycle equations are given by 


0 =1 -Il- 0 1 0. Ojfe] fo 
0 O 1 -1 0 1 Ol}, |=/0 


Multiplying the two matrices on the left, we obtain the equation 


03-6 +04] 10 
=|0 
0 


Ug —Vs +07 
which is equivalent to the three equations on page 74. 


0, +03 +05 +0, 


In this case we wish to eliminate v4, 


So we combine the two equations by 


addition. 
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(b) The chords are AB, CD and CE, and the fundamental cycles are 
ABDEA, BCDB and BCEDB. The fundamental cycle matrix is 


therefore 
BC BD DE EA AB CD CE 


aero) 2 Ya 7 0 6 
cp} 1-1 0 0 0 1 O 
Cel. VY =t =L 10) 0! 10 2 


and the fundamental cycle equations are given by 


0 1 111 0 OJfv,] fo 
1-1 0001 Off, |=|0 
© Si St the 0 


Multiplying the two matrices on the left, we obtain the equation 


Vg t+ V7 +0, +02 
Vy — Ug + U6 = 
03 — U4 — V7 +5 


0 
0 
0 


which is equivalent to the three equations on page 75. 


Solution 2.5 
(a) 
vertices vertices fundamental fundamental 
branch inX my cutset diagram cutset equation 
> ati 4 
Au B 
--| et a 
AB A B,C,D,E (AB, EA) f | i,-i,=0 
D 
BC A,B CDE (BC, BD, EA) i, +i,-i,=0 
cD A, B,C, E D (BD, CD, DE} i, +i, -i, =0 
CE A.B,C,D E (CE, DE, EA} i, +i,-i, =0 
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The cutset (AB, BC, BD} contains the branches AB'and BC. We 
obtain the current law equation for this cutset by combining the cutset 
equations for these two branches, giving 


ist iy -i, = 0. 
(b) 

vertices vertices fundamental fundamental 
branch inX inY cutset diagram cutset equation 
BC AB, DEB C {BC, CD, CE} i, -i; -i, =0 
BD B,C A,D,E (AB, BD, CD, CE} ig tis +ig-i) =0 
DE B,C,D A, E {AB, CE, DE} is +i,-i, =0 
EA B,C, D,E A (AB, EA} i,-=0 


The cutset (AB, BC, BD} contains the branches BC and BD. We 
obtain the current law equation for this cutset by combining the cutset 
equations for these two branches, giving 


i, +i,-1, =0. 


Solution 2.6 


(a) The branches are AB, BC, CD and CE, and the fundamental cutsets 
are {AB, EA}, {BC, BD, EA}, {BD, CD, DE} and (CE, DE, EA}. The 
fundamental cutset matrix is therefore 

AB BC CD CE BD DE EA 
abn) 1 0 0 0 O O +1 
sc} 0 1 0 0 1 O =-1 
co} oO 0 1031-1 0 
CRINO) 10) -0> Bip tO) 1 = 


and the fundamental cutset equations are given by 


(b) 


10000 0 -1]fi,] fo 
01 0).0 1 GO) ailelelo 
0010 1-1 Olfig{“lo 
00010 1 -1}lis| [o 


Multiplying the two matrices on the left, we obtain the equation 
i,-i] [0 
in ip—&|_f0 
igi i (0 
is+i,-i,] [0 
which is equivalent to the four equations on page 76. 


The branches are BC, BD, DE and EA, and the fundamental cutsets 
are {BC, CD, CE}, {AB, BD, CD, CE}, {AB, CE, DE} and {AB, EA}. 
The fundamental cutset matrix is therefore 


BC BD DE EA AB CD CE 
Bf 1 0 0 0 0 ~-1 -1 
papi)oO 100-1 1 1 
pe oO 0 10-10 1 
EA io” OF dt S10) 


0 
and the fundamental cutset equations are given by 


1000 0 -1 -I}fi,] fo 
O40 0-1 2 sel lo 
0010-1 0 1i/l%| Jo 
0001-1 0 offi] {o 

i 

ig 

i; 


Multiplying the two matrices on the left, we obtain the equation 
4-ic-ig) (0 
ig-t, +ig tis |_| 0 
iy byte 
hae 1 


which is equivalent to the four equations on page 77. 


Solution 2.7 


(a) 


(b) 
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AB BC BD CD CE DE EA 
O10: 30) 30" Qe=y 
a 2 100 10) 16 0 

=U GE a SD 60) 

0 Q. =1" 10) 410) 

0 10) 0) OP =te4 


a 
a 
moow>d 


AB BC BD CD CE DE EA 
Ald 0 (0) (6) 10 0) =1 
B,-3}2 1 10 0 oO 0 
Ora 0 =1 10; 2 ok 100 
FO 0 O O-1-1 1 


(c) 


AB BC CD CE 


Solution 2.8 


< 


We first draw the spanning tree and indicate D as the chosen reference 


vertex, 


Tracing the path from each vertex to the reference vertex, we obtain the 


matrix 


Ss 
i t 

eee 
Seesusae 


conc coro 
onoco once 
moco sooo 


1 Ss 
one 
aw oonrn corn 
ocoonro 
coonro - 
3) oo Sonne 
Pushes: | ene: 
<ctrno mano TSS 
‘ = A 1 
2 aon 
226s uf a 
" z ry 
< Tas rl 
i) i) = 


We have 


as required. 


Solution 2.9 


which agrees with the answer given in the solution to Problem 2.6. 


Also 


and so 
0-1-1 0/1 0 0 
CG=|0 0 2 +110 f 0 
1 346) wien 0 


which agrees with the answer given in the solution to Problem 2.4. 


(b) 
CD =[0 12-1 011 0 O}f1 0 0 O}=[0 000 
' 
0 OY ater Toor at OO 10, 01:0 6 
1 EL 2hho! Ora wor a “oP lororore 
ela 
OU Er 16) 
00-1 1 
=-1 -1 0 -1 
This is a special case of the following general result: for an oriented 
graph with n vertices and m edges, the matrix product C,D," = 0. (that 
is, the zero matrix with m—n +1 rows and n - 1 columns). 
Solution 2.10 
4 
, 
3 5 
D 
i 
(a) (b) 


For graph (a), the chords are AD, AC and CD, and the fundamental cycles 
are ADBA, ACBA and BCDB. 
The fundamental cycle equations are 

0, - 05-0, =0 

U3 -% —0,=0 

Uy +, —v; =0. 
For graph (a), the branches are AB, BC and BD, and the fundamental 
cutsets are {AD, AB, AC}, {AC, BC, CD} and {AD, BD, CD}. 
The fundamental cutset equations are 

i, +i, +i, =0 

i, +i, -i, =0 

i, +i; +i, =0. 
For graph (b), the chords are WY, WX and XY, and the fundamental cycles 
are YZWY, ZWXZ and YZXY. 
The fundamental cycle equations are 

v, +0; +0,=0 

U3 + 0% —U, =0 

0, +0; +05=0. 
For graph (b), the branches are YZ, ZW and ZX, and the fundamental 
cutsets are (YZ, WY, XY}, (WY, ZW, WX} and {WX, XY, ZX}. 
The fundamental cutset equations are 

i, -i, -i, =0 

i; -i, -i,=0 

i, -i, +i, =0. 
Graph (b) is the dual of graph (a). See page 20. 
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Notice that the thick edges in graph (a) form a spanning tree; they 
correspond to the thin edges in graph (b); the thin edges in graph (a) 
correspond to the thick edges in graph (b), which form a spanning tree. 


The fundamental cycle equations for graph (a) correspond to the 
fundamental cutset equations for graph (b), and the fundamental cutset 
equations for graph (a) correspond to the fundamental cycle equations for 
graph (b). (Interchange all the vs and is). 


It can be shown that these results hold generally for the dual of a planar 
network. 


Solution 3.1 ' 
This example is similar to Example 3.2. The matrix equation is 


TWO LOCO 6 10 Rey 109 0, wiOm Oy To KO 
' 
Y F000 0 OOM = moon oO) 0) fo. 
OF 0° 9 30) OF 10, 0) NP os Reon Wo) 0 Fo. 
o oO 0 1 0 0 OF Ol 0 OF -K oO 0 Oo 
OF 0) 0 10) 10) Og WO) AGO Re 0) Fo 
007 0 0) 0) Oo oF 0! Oo) eR OO 
' 
OP 10 tO OO 0h on 0 0) 0) on 0) 40 
f 
1 
Solution 3.2 
1 2 3 
-1/ 1 0 
fund) tal cycles: - 
fundamental cutset: Deli} 1 2) 


Combining these matrices with the matrix on page 44, and re-ordering the _ 
variables to agree with the matrices C, and Dy, we get 


0 0 O02 0 Of[y I 
el {0 U2 0 
Bae ET Saree Pa) 2-0 
a ame | 0 i} | 0 
oe ee i} |.0 
a) iz 0 

tT ee Ba 

H xay 


= 
eocoooco 
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Notice that we have made the current source the first component in the 
column vector. It does not matter where the independent sources appear, as 
long as the network voltages and currents appear in the same order in the 
component, cycle and cutset equations. 


Solution 3,3 
1 '-R2 0 0 0 0 0 Oljy 0 
0 '0-RB 0 0 0 O Ofjw| | 0 
0 '0 0 -Rk 0 0 O Olly 0 
0 /0 0 0 -R5 0 0 O}lo5 0 
0 ‘0 0 0 O-R 0 ‘Olly 0 
0 :0 0 0 O O -Ry O}\y, 0 
ee O10) 0), 100. 0y.00))|,24) =e, 
0 0-0 0 i 0 
0 0 0 0 is 0 
1 ig) |.9. 
0 is 0 
0 ig 0 
0 iy 0 
) i; 0 
Solution 3.4 
(a) 
1 0 0 0 0 0:0 oOlfo] [k] [fo 0 0 0 O O}0 olfy 
0 IL 0 0 0 0 {0 O}}m!_j0),J0 0 0 0 0 0;0 1) )d5 
0 0 1 0 O-Ri0 Olly|"|oj"Jo 0 0 0 0 0:0 O]% 
0 0 0 01C 0/0 Of]4] Lo} [0 0 0 O O O}1 Oh 
iy L 
it ig 
V2 % 
ig iy 


Note that the variables whose derivatives occur in the component 
equations (that is, v, and i,) appear last in the column vectors. 


(b) The fundamental cycle and cutset matrices can be rewritten as 


follows: 
1100 0 0}1 O7f?, 0 Note that we have rearranged the 
0 10 ! columns of the matrix so as to put the 
=t 0 0/0 0 ||?) _10 edge variables in the same order as in 

0001-1010 O})% 0 the component equation matrix. 
0000-1110 1}/4]} [o 

h 

i, 

2% 

Cy 


Combining this with the component-equation matrix, we obtain the 
H-matrix equation Hx = y + Kx, as follows: 


1 0 0 0 0 0:0 Offm] [kK] fo 0 0 0 0 O}0 Offe, 
0 WL 0 0 0 0:0 f/m] jo} jo 0 0 0 0 0:0 1I/s, 
0 0 1 0 O-R4i 0 O]]m) Jo} |O 0 0 0 0 0:0 Oly 
0 0 0 0 YC 0:0 Of]4| Jo) JO 0 0 O O O:1 OA 
1 1 0 0 O O12 Ol#lelOj,Jo 0 0 0 0 O10 Oh 
0-1 1:0 0 0:0 O/m] jo} jo 0 0 0 O 0:0 O]% 
0 0 0 1-1 0:0 Oj}m] Jo} jo 0 0 0 0 0/30 Offe, 
0 0 0 0-1 1:0 lis} lol] [oO 0 0 0 O O'0 Oh 
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(c) To find the state equations we start with rows 5 and 8 and perform 
Tow operations in such a way that the first six entries in these rows 
are all zero. 


First state equation: Take row 5, and subtract row 1 (to eliminate the 
1in column 1) and L x row 2 (to eliminate the 1 in column 2), 
Performing the same operations on the right-hand side of the 
equation yields the state equation 

2, =-k-Li; 


Second state equation: We can eliminate the entries in columns 5 and 
6 of row 8 by taking 


(row 8) + (C x row 4) + ((1/R,) x row 3)—((1/R,) x row 6) 
—((L/R,) x row 2) 


This yields the state equation 


: «Ly 
I= CO; Rta? 


Index 
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3-terminal 6 
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controlled source 14 
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current source 6 

cutset rank 31 
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dual network 20 
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electrical network 6 
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fundamental cutset 30 
equation 31 

matrix 32 
fundamental cycle 25 
equations 27 

matrix 28 


Gausssian elimination 50 
graphical representation 15 

of 2-terminal component 10 

of multi-terminal component 15 


Hooke’s law 9, 13, 18 
hydraulic network 7, 9 


ideal independent source 13 
ideal transformer 14 
incidence matrix 
of oriented graph 34 
reduced 34 
independent source 13 
inductance 12 
inductor 12 
inverse matrix 35 


Kirchhoff’s current law 8, 15 
Kirchhoff's voltage law 8, 15 


linearity 22 


mechanical system 7 
multi-terminal component 10, 14 


network 

dual 20 

electrical 7 
Newton’s second law 13 
n-terminal component 17 


Ohm’s law 9, 12 
oriented graph 5, 11 
reduced incidence matrix 34 


partitioned matrix 35 
potential difference 7 


potential difference conservation law 8 


prescribed current 13 
prescribed voltage 13 
pressure difference 9 
principle of superposition 22 


reduced incidence matrix 34 
reference direction 10 
reference vertex 34 
resistance 12 

resistor 6, 9, 10, 12 


spanning tree 26 

sparse matrix 45 

spring 7,9 

state equation 42, 55 

state variable 42, 54 

system independent behaviour 22 


Tellegen’s theorem 38 
terminal 6 

2-terminal component 6, 10 
3-terminal component 6 
4-terminal component 6 
through variable 5, 7 
time-invariance 22 
transformer 6, 14 

transistor 6, 14 


variable 
across 7 
state 42, 54 
through 7 
vertex law 15 
equation 9 
voltage 7 
voltage source 6 
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